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SUMMARY 

This  document  contains  a  review  of  the  Workshop  proceedings,  together  with  the  extended 
abstracts  submitted  by  the  speakers. 
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Finite  Elements  in  Composite  Materials  Workshop 

Papers 


Paper 

1  Prof.  F-K  Chang,  Stanford  University 

Analysis  of  laminated  composites  beyond  initial  failure. 

2  Mr.  J.  Underwood,  US  Army  R&D  Centre. 

FE  calculation  of  stresses  and  displacements  in  composite  cylinders. 

3  Prof.  O.  Ochoa,  Texas  A&M  University 
Failure  mechanics  modelling  using  ABAQUS. 

4  Dr.  M.  Wisnom,  University  of  Bristol 

Using  ABAQUS  to  model  composite  structures  and  materials. 

5  Dr.  W.S.  Arnold,  Paisley  College 

Some  considerations  on  the  use  of  FE  with  composite  structures. 

6  Mr.  M.J.  Gunn  and  Prof.  L.  Holloway,  University  of  Surrey 
Finite  element  analysis  of  composite  space  systems. 

7  Mr.  J.D.  MeVee  and  Mr.  R.S.  Dow,  ARE,  Dunfermline 
Finite  clement  analysis  of  composite  naval  structures. 

8  Mr.  P.  Dallard,  Ove  Amp 
Analysis  of  composite  racing  yachts. 

9  Prof.  G.  Verchery,  Ecole  des  Mines,  St.Etienne 

Some  developments  in  finite  element  aided  design  for  composite 
stmetures. 

10  ♦  Dr.  A.  Miravete,  University  of  Zaragoza 

FEM  applied  to  a  tmek  frame  made  of  hybrid  composites. 

11  Dr.  H.  Schellekens,  Technology  University,  Delft 

Numerical  problems  encountered  in  non-linear  finite  element  analysis 
of  free  edge,  delamination. 

12  Dr.  R.  Rolfes,  DIT.R,  Braunschweig 

Thermal  analysis  of  laminated  composites;  element  development 
and  experience  with  MSC/NASTRAN. 

1 3  Mr.  R.  Scherer,  University  of  Hamburg-Harburg 

Simulation  of  thermoforming  process  of  continuous  fibre-reinforced 
thermoplastic  using  MARC. 

*  Dr.  Miravete  was  unable  to  present  his  paper.  An  additional  paper  was 
read  in  its  place  by  Dr.  N.  l^lan,  ICl. 
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INTRODUCTION 


1 . 


The  objective  of  this  Workshop  was  to  consider  the  current  position  of  applying  the  finite 
element  method  to  applications  involving  composite  materials.  It  was  expected  that  issues 
needing  attention  would  be  identified.  ^  / .  •  ^  •.  z  ^ 

To  achieve  these  ends  13  speakers  representing  several  aspects  of  the  field,  from  a 
number  of  countries  (France,  1 :  Germany,  2;  Netherlands,!;  Spain,!;  UK,  5;  USA,  3)  were 
invited.  All  accepted,  but  one  withdrew  just  before  the  meeting.  This  talk  was  replaced  by  a 
contribution  from  one  of  the  delegates.  There  was  one  further  contributed  paper. 

Seven  software  companies  initially  agreed  to  participate,  via  discussions  and 
demonstrations;  finally  only  three  attended.  Some  withdrew  because  they  wanted  significantly 
more  time  for  demonstrating  than  was  scheduled  in  the  programme. 

In  total  52  persons  attended  the  workshop  made  up  as  follows;  !3  speakers,  28 
delegates,  3  representatives  from  software  suppliers,  8  Imperial  College  staff  (see  Appendix 
C). 

The  meeting  was  chaired  by  Dr.  I.  C.  Taig,  ex-British  Aerospace,  Consultant  to  NAFEMS 
(National  Agency  for  Finite  Element  Methods  &  Standards). 

PROCEEDINGS 

The  format  of  the  workshop  allowed  for  10  minules  discussion  after  each  presentation, 
with  additional  discussion  periods  scheduled  at  the  end  of  the  first  day  and  at  the  end  of  the 
meeting.  To  keep  within  time  constraints  discussion  always  had  to  be  stopped  after  each  paper, 
although  the  general  discussion  periods  appeared  to  be  of  adequate  length  to  allow  all  present  to 
make  their  points. 

Extended  abstracts  of  the  invited  papers  and  original  contributed  paper  are  included  as 
Appendix  A.  The  substituted  paper  (given  subsequently  at  the  ICCS-6  Conference)  is  included  as 
Appendix  B. 

In  his  opening  address  Dr.  Taig  noted  that  none  of  the  currently  available  finite  element 
systems  gave  a  true  representation  of  a  comfxtsite.  The  requirements  could  be  split  into  four 
categories  (essential,  expected,  desirable,  special)  which  can  then  be  addressed  under  a  number 
of  areas:  materials,  geometry,  formulation,  solution,  results,  failure  [Ref:  'Finite  Element 
Analysis  of  Composite  Materials',  I.  C.  Taig,  NAFEMS,  East  Kilbride,  Scotland,  1990]. 

Many  of  the  speakers  addressed  most  of  the  points  listed  by  Or.  Taig.  In  particular  the 
following  were  emphasised:  the  need  for  an  accurate  description  of  the  material;  the  need  to 
model  the  failure  process;  the  need  for  systems  that  are  easy  to  use;  the  fact  that  different 
packages  can  produce  different  results. 

The  conclusions  from  the  discussions  and  presentations  are  given  in  the  next  section. 
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CONCLUSIONS 


The  workshop  was  generally  agreed  to  have  been  a  success  and  to  be  worth  repeating. 

The‘‘technical  conclusions  were  as  follows; 

Materials  "<0  -ttie  actual  material  must  be  modelled - 

■  <  > properties  must  be  appropriate  to  the  scale  of  the  analysis: 
Geometry  -■  -  the  system  must  allow  for  correct  orientation  of  the  material; 
independent  of  the  mesh 
it  should  be  possible  to  model  realistic  lay-ups 
-  it  should  be  possible  to  represent  (ply)  discontinuities 

Formulation  it  should  be  possible  to  incorporate  through-thickness  modelling. 

^  s  specific  formulation  is  needed  to  allow  for  iteration : 

Solution  -  stable  methods  are  needed  to  model  damage  propagation 

^  ^  it  should  be  possible  to  include  material  and  geometric  non-linearity 

Results  -  need  transverse  and  through-thickness  (as  well  as  In-plane)  results 

need  flexible  post-processing,  with  ability  to  discard  irrelevant 
information 

Failure  ’  include  relevant  critical-  properties ; 

r  system  must  be  able  to  identify  failure  mode(s);  ■, 

must  be  able  to  model  progression  of  failure  to  final  collapse 
General  '■  consider  using  spreadsheets  to  interface  input  data  with  FE  system 
users  need  more,  and  better,  communication  with  software  vendors 
establish  a  system  (benchmarks)  for  verifying  results  .  .. ,, 


RECOMMENDATIONS 

the  workshop  organisers  to  consider  repeating  the  event  in,  say,  two  years'  time, 
attendees  to  be  sent  a  questionnaire  by  the  organisers  to  seek  views  on  fomnat  and 
content  of  a  future  meeting. 

software  vendors  and  NAFEMS  to  be  sent  a  copy  of  the  current  document  in  order 
to  encourage  action  on  the  conclusions  of  the  workshop. 


APPENDIX  A 


We  are  all  familiar  with  the  potential  of  both  composite 
materials  and  finite  element  analysis.  It  is  difficult  to  realise  the 
full  potential  of  a  composite,  and  there  must  be  many  instances 
of  such  materials  being  abandoned  in  favour  of  metals 
because  of  a  lack  of  understanding  of  some  fundamental 
aspect  of  the  composite’s  performance.  Equally,  incorrect 
answ/ers  with  consequent  reflection  on  the  behaviour  of  the 
product,  can  be  obtained  from  a  finite  element  analysis  if 
insufficient  attention  is  paid  to  details  of  the  modelling. 

When  it  comes  to  using  finite  elements  with  composites  we 
have  at  our  disposal  greater  potential  but  more  scope  for 
making  mistakes.  In  view  of  the  increasing  use  of  both  areas 
of  technology  it  seemed  timely,  therefore,  to  hold  a  workshop  to 
discuss  the  issues  involved,  to  identify  problem  areas,  and  to 
recommend  ways  of  improving  our  use  of  F-E  with  these 
materials. 

In  putting  this  meeting  together  we  tried  to  assemble  speakers 
who  are  practitioners,  and  can  thus  talk  about  real  problems 
and  how  they  solved  them,  in  several  fields  of  application, 
together  with  software  suppliers.  The  workshop  will  only 
succeed  if  we  have  a  dialogue  between  users  and  suppliers. 

We  are  extremely  grateful  for  the  financial  support  provided  by 
the  European  Research  Offices  of  the  United  States  Army  &  Air 
Force.  We  are  particularly  grateful  for  the  encouragement 
provided  in  the  early  planning  stage  by  Dr.  Wilbur  Simmons  of 
the  US  Army  Office. 

F.L.  Mathews, 

Director,  Imperial  College  Centre  for  Composite  Materials. 
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ANALYSIS  OF  LAMINATED  COMPOSITES  BEYOND  INITIAL  FAILURE 


WORKSHOP:  FINITE  ELEMENT  METHODS  FOR  COMPOSITES 

Imperial  College  of  Science  Technology  and  Medicine 
University  of  London 
September  5-6,  1991 


Fu-Kuo  Chang 

Department  of  Aeronautics  and  Astronautics 
Stanford  University,  Stanford,  CA  94305 
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ABSTRACT 


In  order  to  analyze  the  mechanical  response  of  fiber-reiiiforccd  laminated  composites,  an 
analytical  model  is  required  which  can  calculate  accurately  stresses,  strains  and  deformations 
of  the  composites,  and  predict  reliably  the  state  of  damage  and  the  materia!  p'opcrties  at 
every  damaged  state  inside  the  materials  as  a  function  of  the  applied  load  and/or  the  history 
of  the  load.  Due  to  inhomogenity  and  anisolrophy,  stress  distributions  inside  composites  are 
very  complicated  and  may  vary  from  layer  to  layer  and  location  to  location.  As  a  c  snsequence, 
damage  in  composites  involves  multiple  failure  modes,  which  are  strongly  influen  ;cd  by  stress 
and  strain  distributions  inside  the  materials. 

Finite  element  methods  have  been  developed  and  utilized  in  the  literature  for  calculating 
stresses  and  deformations  of  laminated  composites  before  the  occurrence  of  damage.  Once 
the  materials  suffer  some  degree  of  damage,  an  appropriate  failure  analysis  is  needed  in 
conjunction  with  the  finite  element  analysis  for  analyzing  the  post-damage  response  of  the 
composites.  In  general,  failure  modes  in  laminated  composites  can  be  classified  into  two 
types:  in-plane  failure  such  as  fiber  breakage,  fiber  compresssion,  matrix  cracking,  matrix 
compression,  fiber  matrix  shearing,  and  out-of- plane  failure  such  as  delamination.  The  failure 
mechanism  of  each  failure  mode  is  different;  hence,  the  effect  of  damage  on  the  material 
properties  is  also  strongly  dependent  upon  the  mode  of  failure. 

Several  methods  have  been  proposed  in  the  literature  for  predicting  failure  of  laminated 
composites  beyond  initial  failure.  In  this  presentation,  a  progressive  failure  analysis  based 
on  the  continuum  mechanics  concept  will  be  presented  for  predicting  the  in-plane  damage 
in  composites.  The  method  consists  of  a  nonlinear  finite  element  analysis  for  calculating 
stresses,  and  a  failure  analysis  for  predicting  the  type  and  the  extent  of  damage  in  the  mate¬ 
rials  as  a  function  of  the  applied  load.  The  application  of  this  analysis  will  be  demonstrated 
through  modeling  the  response  of  laminated  composites  containing  holes  and  cutouts. 

For  analyzing  delamination  failure,  a  model  will  be  presented  which  consists  of  a  non¬ 
linear  finite  eleinent  method  for  calculating  stresses,  a  contact/slip  analysis  for  modeling  the 
interfacial  condition  of  the  delamination  surfaces  during  loading,  and  a  crack  growth  criterion 
for  predicting  delamination  propagation.  The  model  will  be  demonstrated  for  modeling  the 
delamination  propagation  in  composites  resulting  from  compression  or  transverse  loading. 
Discussion  will  also  be  given  on  the  limitation  of  the  application  of  both  analyses  and  the 
current  efforts  for  improving  the  analyses. 
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FINITE  ELEMENT  RESULTS  FOR  COMPOSITES  IN  ARMRMENT  APPLICATIONS! 
PRESSURIZED  CYLINDERS  AND  FRACTURE  TEST  METHODS 


by  J.H.  Underwood  and  M.D.  Witherell 

US  Army  Armament  Research,  Development  and  Engineering  Center 
Benet  Laboratories  Watervliet,  LY  12189  USA 


Hollow  cylinder  applications  far  composites  will  be  discussed, 
as  related  to  armament  components.  Results  oT  recent  and  ongoing 
stress  analysis  af  various  mul  t )  or thotr op i c -1 ayer ed  cylinders  will  be 
described,  simulating  the  type  of  applied  and  residual  loads  common 
to  thick-wall  cylinders  and  the  material  properties  of  both  metallic 
and  composite  materials. 

The  need  for  fracture  mechanics  analysis  of  damaged  composite 
structures  will  be  proposed,  particularly  in  relation  to  armament. 

The  various  ASTM  efforts  to  develop  fracture  tests  of  composites  will 
be  summarized,  with  emphasis  on  finite  element  calculations  of  stress 
intensity  factors  for  laminates  in  various  test  configurations. 


for  presentation  at: 

Workshop  on  Finite  Element  Methods  for  Composites 
Imperial  College  London  5/6  Sept  1991 
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F INITE  ELEMENT  RESULTS  FOR  COMPOSITES  IN  ARMAMENT 
APPLICATIONS;  CYLINDERS  AND  FRACTURE  TOUGHNESS 

J.H.  Underwood  and  M.D.  Witherell 
Army  Armament  RD&E  Center  Watervliet,  NY  12189 

COMPOSITE  JACKETTED  CYLINDERS: 

■  aPDlicatlon:  light  weight  Jacket  over  alloy  steel 

-  titanium  alloy;  brittle  I 

-  carbon/polymer;  higher  atiffneaa/weight  and  toughness 

•  multi-laver  composite  cylinders 

-  use  a  modified  Lekhnitskii  analysis 

-  close  agreement  with  ABACUS  results 

-  can  model  Interference  residual  stress 


TRANSLAMINAR  FRACTURE  TOUGHNESS: 

■  aoDllcatlon:  composite  structures  on  battlefield 

~  projectile  impact  /  penetration 

-  translaminar  fracture  vs  interlaminar  delamination 

■  Harris  A  Morris  approach 

-  FE  calculations  of  K-applied  at  notch 

'  tension,  compact  and  bend  specimen  confifguratlons 

-  orthotropic  model  of  carbon/polymer  laminates 

-  K-ortho  •  K-iso  ;  except  In  tension 


5 


Hfitr  BTlilV  (lilCIO  IK. /IK.) 


RATIO  OF  STRESS-INTENSITY  FACTORS  FOR  ORTHOTROPIC  AND 
ISOTROPIC  MATERIAL  FROM  FINITE-ELEMENT  ANALYSIS 


1.2 


K 


ortho 


1.0 

0.8 

0.6 

0.4 

0.2 


^  9  8  Q 
□  □  □  □ 


Middle— crock  tension: 
O  [0/±45]^3 
a  [0/90]ns 
A  [0/±45/90]„3 


0.0  i - ^ 

0.0  0.2 


0,4  0.6 

o/w 


0.8 


1.0 


DELAMINATION  FRACTURE  MODES 


Mode  1  Mode  II  Mode  HI 


interlaminar 

tension 
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Translaminar  Fracture  /  E24. 07.02 

PLANS:  Specimen  Configurtations 
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FAILURE  MECHANISMS  AND  MODELLING 


Ozden  O.  Ochoa 

Center  for  Mechanics  of  Composites 
Texas  A  &  M  University 
College  Station,  Texas  77843-3123 

ABSTRACT 

Structural  response  of  composite  components  presents  a  major  challenge  in 
assessing  and  modelling  damage  mechanisms  as  a  function  of  load.  Quite  often 
damage  modes  are  present  as  a  byproduct  of  manufacturing  cycles  prior  to  any 
structural  loading.  The  task  of  identifying  the  present  anomalies  and  numerically 
capturing  their  growth  is  of  paramount  importance  and  difficulty. 

Following  is  a  brief  discussion  on  three  different  problems  to  illustrate  a 
practical  procedure  in  integrating  present  defects  in  numerical  approximations. 
These  cases  represent  modelling  of  a  cutout  section  in  the  web  of  a  stiffened  panel, 
a  loaded  hole  with  drilling  defects  and  a  unit  cell  concept  associated  with  textile 
architecture  of  woven  prepreg. 


Cutouts  in  a  Stiffened  Panel 

An  area  of  critical  concern  for  load-carrying  members  is  cutouts  which  may 
provide  a  passage  for  hydraulic  lines  or  avionic  harnesses,  or  on  a  larger  scale,  an 
access  door  in  an  aircraft  fuselage.  Parametric  studies  are  undertaken  to  reflect  the 
differences  in  the  selection  of  stacking  sequences,  cutout  size,  and  geometry.  Primary 
emphasis  is  on  the  identification  of  localized  high  stress  and  strain  in  order  to  warn 
the  designer  of  potential  hazards.  In  this  study,  Patran  II  software  is  used  to  develop 
the  nodal  and  element  geometry,  then  converted  to  a  neutral  file.  The  neutral  file 
is  used  as  input  to  Abaqus,  a  general  purpose  code  that  is  a  robust  performer  in 
nonlinear  behavior. 

The  response  of  a  hat  stiffened  panel  with  a  cutout  in  the  web  section 
subjected  to  bending  loads  is  considered.  Since  tensile  loads  are  primarily  sustained 
by  the  flat  panel  and  not  transferred  into  the  stiffener,  out  of  plane  loading  is 
selected  to  observe  the  effect  of  loads  on  the  web.  Therefore,  bending  moments 
along  the  panel  axis  are  imposed  to  induce  a  stress  field  across  the  cutout  region. 
This  loading  is  representative  of  the  in-flight  loading  incurred  by  a  stiffened  wing 
skin.  Shell  elements  are  used  in  the  modelling  throughout  the  structure.  The  skin 
configuration  is  a  [-45/45 /90/-45/0/45/0/-45/0/45/0]s  laminate.  The  stiffener  web 
is  three  layers  of  ±45®  fabric,  with  a  12.7  mm  diameter  cutout  through  each  web 
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center.  The  cutout  is  positioned  79.4  mm  from  each  panel  edge,  and  is  located  25.4 
mm  vertical  from  the  panel.  The  flange  is  reinforced  with  0°  plies,  i.e.;  (45/- 
45/04/45/-45/0^4S/-45]. 

The  web  section  predominantly  carries  shear  stresses  throughout.  The  stress 
results  for  circular  and  a  diamond  cutout  in  the  web  supports  the  use  of  diamond 
shaped  cutouts  in  regions  of  a  structure  which  are  primarily  in  shear  loading,  and 
contain  primarily  ±45°  fibers.  It  can  be  deduced  that  for  an  effective  composite 
design,  a  dominant  percentage  of  fibers  should  be  aligned  with  the  principle  stress 
direction.  In  addition,  if  a  cutout  is  required  in  the  design,  it  is  most  desirable  to 
implement  a  rectangular  cutout  which  contains  sides  parallel  and  perpendiai'ar  to 
the  principle  stress  directions. 


Hole  Quality  and  Perfonnance  in  Joints 

Frequently,  drilling  causes  delaminated  regions  around  the  circumference  of 
the  hole.  The  propagation  of  these  interlaminar  cracks,  or  delaminations,  is  one  of 
the  most  serious  problems  in  failure  prediction  of  composite  parts.  Actual  damage 
zones  caused  by  drilling  can  be  documented  through  the  use  of  both  X-ray  and 
ultrasonic  C-scans  and  optical  micrographs.  From  the  study  of  optical  micrographs, 
it  is  observed  that  the  lower  30%  to  40%  of  the  coupon,  drill  exit  side,  contains 
delaminated  layers.  Therefore  in  the  numerical  model,  they  will  be  assigned  reduced 
moduli.  Furthermore,  as  an  approximation,  we  will  treat  the  damage  regions  as 
circular  in  shape  and  centered  about  the  center  line  of  the  hole  with  effective 
delamination  diameters  of  1.5-D  and  2.0-D. 

A  pin  loaded  coupon  under  tension  load  is  simulated  with  a  finite  element 
model.  TTie  pin  is  represented  by  distributed  loads  applied  to  its  contact  surface. 
Thus  rather  than  model  the  pin  itself,  the  effects  of  the  pin  are  represented  as  an 
evenly  distributed  pressure.  Since  the  coupon  geometry,  assumed  load  distribution, 
and  assumed  delamination  shape  are  symmetrical  about  the  long  axis  of  the 
specimen,  a  half  symmetry  model  is  generated.  Finite  element  analysis  is  performed 
microcomputer  with  a  MS  DOS-based  program.  Algor.  A  four  node  composite  plate 
element  based  on  the  assumptions  of  classical  plate  theory  is  used.  Each  node  has 
five  degrees  of  freedom:  two  orthogonal  in-plane  displacements  (u,v);  one  out-of¬ 
plane  displacement  (w);  and  rotations  about  two  orthogonal  in-plane  axes  (^*w/^x^ 
dSv/5y^). 

Despite  some  short  conungs  in  the  modelling  assumptions,  these  two- 
dimensional  models  give  insight  into  the  conditions  surrounding  the  hole.  Through 
a  layer  by  layer  representation  of  the  lamina  stress,  these  models  give  a  better 
understanding  of  the  damage  created  in  the  coupons  during  testing.  The  results 
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describe  the  stress  distribution  around  the  hole  allowing  comparisons  between 
stacking  sequence,  load  levels  and  drilling  induced  damage  levels. 


IVoven  Prepreg/Laminate 

In  brittle  material  systems,  parameters  of  importance  in  modelling  constitutive 
behavior  are  largely  driven  by  processing  methods.  For  example,  shrinkage  observed 
during  processing  results  in  a  laminate  vhich  exhibits  both  matrix  cracks  and  weak 
interfacial  bonds.  It  is  important  to  assess  the  status  of  a  laminate  prior  to 
mechanical  and  thermal  service  loads.  However,  it  may  be  improper  to  assume  that 
the  above  cited  damage  is  always  detririental.  A  strong  emphasis  need  to  be  placed 
on  obtaining  statistical  data  to  support  the  modelling  of  possible  defects  within  a  unit 
cell  of  the  woven  laminate.  The  goal  is  to  illustrate  the  changes  in  the  moduli 
(constitutive  relations)  as  a  function  of  applied  load. 

The  numerical  models  that  follow  address  the  response  of  thin  laminates,  i.e., 
six  layers  of  8-Harness  satin  weave  woven  laminates,  subjected  to  static  loads.  After 
careftil  studies  of  the  surface  with  SEM,  the  unit  cell  is  approximated  as  a  simple 
square  with  2300  micron  edges.  A  bundle  crossover  square  region  of  800  microns  in 
the  center  is  modelled  with  a  typical  adjacent  debond  of  length  240  microns.  Since 
the  material  in  the  immediate  vicinity  of  the  crossover  region  consist  of  warp  and  fill 
yarns,  it  is  treated  as  a  homogeneous  specially  orthotropic  material.  The  material 
properties  of  bundle  crossover  section  is  also  treated  as  specially  orthotropic  , 
however  they  are  modified  as  a  function  of  the  crimp  angle. 

Four  node  shell  elements,  S4R,  of  ABAQUS  with  nonlinear  geometric  option 
and  a  user  material  subroutine  developed  by  the  author  is  used.  This  specially 
developed  subroutine  enables  the  incremental  analysis  to  check  for  failure  strain  in 
each  element,  at  each  integration  point  and  at  each  load  level.  For  the  tensile  loads, 
symmetry  boundary  conditions  are  applied  to  a  quarter  model.  The  tensile  stress- 
strain  results  indicate  that  the  load  carrying  ability  is  significantly  diminished  in  the 
presence  of  a  debond.  Since  the  actual  laminate  consists  of  many  layers  of  the 
woven  fabric  with  random  defects,  it  is  appropriate  to  calculate  an  effective  elastic 
modulus  from  the  unit  cell  analysis  that  can  be  used  as  a  layer  modulus  in  structural 
response  models. 
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USE  OF  ABAQUS  TO  MODEL  COMPOSITE  STRUCTURES  AND  MATERIALS 


Michael  R.  Misnom 

University  of  Bristol,  Department  of  Aerospace  Engineering 


1 .  Introduction 

ABAQUS  has  been  used  at  Bristol  University  for  a  wide  range  of 
composite  problems  from  analysis  of  complete  structures  down  to 
micromechanical  modelling  at  the  level  of  individual  fibres.  In 
this  paper  a  number  of  applications  are  presented  and  the  way  that 
ABAQUS  was  used  is  discussed.  Special  features  that  make  AEAQUS 
suitable  for  the  analysis  of  composites  are  described,  as  well  as 
some  of  the  limitations  and  problems  encountered. 


2.  Cruciform  Joint 

The  first  application  is  the  analysis  of  thermoplastic 
matrix-carbon  fibre  cruciform  joints.  This  was  a  conventional 
structural  analysis  which  illustrates  some  of  the  features  ABAQUS 
has  for  modelling  composites.  Layered  shells  were  used,  and  these 
are  available  as  4  noded  or  8  noded  elements.  Many  finite  element 
programs  calculate  equivalent  orthotropic  properties  for  membrane, 
bending  and  membrane-bending  coupling,  and  this  limits  them  to 
linear  elastic  analysis.  However  ABAQUS  uses  numerical 
integration  through  the  thickness,  which  allows  non-linear 
material  response  to  be  modelled  in  the  layers.  These  elements 
also  allow  transverse  shear  flexibility. 

Linear  elasticity  was  assumed  in  this  particular  analysis.  A 
static  analysis  was  carried  out  of  the  joint  under  tension 
loading.  Thermal  residual  stresses  were  included  by  imposing  a 
temperature  change.  The  analysis  was  relatively  straight  forward, 
and  no  particular  problems  were  encountered.  A  large  amount  of 
output  data  was  produced  of  layer  stresses  and  strains.  These  can 
be  output  at  integration  points  or  averaged  at  nodes.  One 
disadvantage  is  that  Interlaminar  shear  stresses  are  not  computed. 
However  shear  forces  can  be  recovered  in  the  output  file  and 
interlaminar  stresses  calculated  in  a  post-processing  routine. 
Another  limitation  is  that  standard  failure  criteria  are  not 
available,  although  again  these  can  be  included  in  the  post 
processing. 


3.  Pin  Loaded  Plates 

Another  application  is  in  the  analysis  of  failure  of  pin-loaded 
silicon  carbide-aluminium  laminated  plates  loaded  in  tension  [1]. 
The  local  stress  distribution  is  determined  by  the  contact  between 
the  pin  and  the  hole.  ABAQUS  has  conventional  gap  elements  to 
model  contact,  however  it  also  has  interface  elements  which  are 
particularly  useful.  The  advantage  of  these  elements  is  that  it 
is  not  necessary  to  use  local  coordinate  systems  to  define  the 
direction  of  the  interface.  Also  interface  stresses  are 
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calculated  rather  than  forces 


The  stress-strain  response  of  the  material  Is  non-llnear, 
especially  In  shear.  This  irfas  modelled  using  the  user  defined 
subroutine  UMAT.  For  a  given  strain  vector  it  is  necessary  to 
define  the  tangent  stress-strain  matrix  and  the  stress  vector. 
The  stress-strain  response  of  the  unidirectional  material  in  the 
transverse  and  fibre  direction  in  both  tension  and  compression  and 
also  the  shear  response  were  measured  and  the  data  was  curve 
fitted.  The  orthotropic  stress-strain  matrix  was  then  calculated 
by  interpolating  the  data.  The  response  was  assumed  to  be 
elastic.  In  practice  there  is  plastic  deformation,  but  for 
nonotonically  Increasing  load  this  will  not  affect  the  results. 
This  approach  assumes  that  the  response  in  the  transverse 
direction,  fibre  direction  and  in  shear  are  independent.  For 
example  the  non-linearity  in  transverse  compression  is  assumed  not 
to  be  affected  by  the  shear  stress.  In  practice  there  is  some 
interaction,  but  this  is  a  simple  approach  w.iich  allows  a 
reasonable  approximation  of  the  complex  behaviour  of  the  material. 

Within  the  UMAT  routine  it  was  also  possible  to  include  a  number 
of  different  failure  criteria.  Partial  failure  could  also  be 
modelled.  For  example  when  the  transverse  tensile  stress  in  one 
layer  reached  a  critical  value  the  transverse  an'^  shear  moduli 
could  be  set  to  zero,  whilst  still  retaining  the  iit/re  direction 
modulus.  This  allowed  failure  to  be  tracked  from  initiation  right 
through  until  catastrophic  fracture.  Comparison  with  test  results 
indicated  a  reasonable  correlation. 


4.  Four  Point  Bending  Tests 

Four  point  bending  of  unidirectional  carbon  fibre-epoxy  was 
analysed  in  order  to  quantify  the  errors  involved  in  applying 
linear  elastic  bending  theory  to  deduce  the  stress  at  failure  from 
the  applied  load  [2].  A  two  dimensional  analysis  of  a  slice 
through  the  specimen  was  carried  out.  Very  large  deformations 
occur  in  these  tests  and  so  the  non-linear  geometry  option  was 
used.  Due  to  the  large  rotations,  the  angle  of  the  forces  at  the 
rollers  changes  significantly.  Also  the  specimen  slips  on  the 
rollers  and  so  the  point  of  contact  moves  along  the  specimen. 
Both  of  these  effects  can  be  modelled  using  the  rigid  surface 
capability  in  ABAQUS.  The  rollers  are  defined  as  rigid  surfaces 
and  interface  elements  are  used  to  connect  them  to  points  on  the 
specimen  where  contact  may  occur.  At  each  iteration  all  these 
Interfaces  are  checked  for  contact  and  so  the  actual  points  of 
contact  do  not  have  to  be  defined  beforehand.  This  worked  very 
well,  although  it  was  found  that  a  very  fine  mesh  was  needed  in 
order  to  avoid  step  changes  in  results  as  contact  shifted  from 
node  to  node  along  the  specimen. 

Friction  was  also  Included  at  the  rollers.  This  worked  well  for 
small  values  of  friction  coefficient,  however  for  larger  values 
convergence  became  slow.  The  friction  algorithm  makes  use  of  a 
stiffness  in  stick,  the  value  of  which  is  hard  to  determine.  If 
too  low  a  value  is  used  significant  relative  movement  can  occur 
whilst  the  interface  is  supposed  to  be  sticking.  Alternatively, 
if  too  high  a  value  is  used  convergence  can  become  difficult. 
This  problem  should  be  solved  shortly  with  the  introduction  of  a 
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new  friction  algorithm. 

Non-linear  material  properties  were  required,  in  this  case  to  take 
account  of  the  changing  modulus  with  strain  in  the  fibre 
direction.  This  was  modelled  with  an  orthotropic  non-linear 
elastic  UMAT  routine,  similarly  to  the  metal  matrix  plate 
analysis.  Because  of  the  large  displacements,  rotation  of  the 
material  axis  system  with  the  specimens  is  essential.  ABAQUS 
handles  this  provided  the  ^ORIENTATION  option  is  used.  Stresses 
are  also  output  in  this  rotated  system.  Strains  of  over  2%  arise 
in  tests  on  carbon  fibre-epoxy  and  with  glass  fibre-epoxy  strains 
of  over  5%  can  occur.  Whilst  ABAQUS  can  handle  such  large 
strains,  it  should  be  noted  that  errors  can  arise  in  using  a 
standard  elasticity  formulation  at  large  strains.  These  errors 
are  typically  of  the  same  order  of  magnitude  as  the  strain. 
Continuum  elements  in  ABAQUS  use  true  strain  rather  than 
engineering  strain  and  at  large  strains  the  difference  between 
these  strain  measures  also  starts  to  becone  significant. 

This  analysis  generally  worked  very  well.  However  one  problem 
arose  when  a  similar  analysis  was  performed  with  shell  elements. 
It  was  found  that  under  pure  bending  small  fictitious  axial  forces 
were  generated.  The  same  phenomenon  occurred  with  beam  elements. 
Although  the  magnitude  of  the  forces  was  small,  they  were 
significant  since  the  whole  purpose  of  the  analysis  was  to 
investigate  deviations  from  standard  bending  theory.  This 
illustrates  the  care  that  has  to  be  taken  in  any  analysis  to 
adequately  check  and  understand  the  procedures  being  used  by  using 
simple  test  cases. 


5.  Delamination  of  Tapered  Specimens 

Oelamination  in  tapered  unidirectional  glass  fibre-epoxy  with 
dropped  plies  has  been  analysed  [3}.  A  plane  stress  analysis  of  a 
slice  through  the  specimen  was  carried  out  using  orthotropic 
elastic  properties.  In  the  taper,  the  fibres  curve  around  the  end 
of  the  terminating  plies  and  so  there  is  a  continuous  change  of 
fibre  direction.  This  effect  can  be  modelled  very  easily  using 
the  user  subroutine  ORIENT  to  define  the  material  axes,  but  would 
cause  difficulties  without  this  option. 

The  stresses  in  the  specimen  subject  to  tension  were  calculated 
and  then  the  strain  energy  release  rate  was  evaluated  assuming 
delamination  above  and  below  the  terminating  plies  propagating 
into  the  thick  section.  The  energy  release  rate  was  calculated 
from  the  difference  in  the  total  strain  energy  in  analyses  before 
and  after  propagation  of  the  delamination. 

It  has  been  found  that  the  resin  layer  between  plies  affects  the 
delamination  significantly.  This  effect  has  been  modelled  using 
non-linear  springs  to  represent  the  elastic-plastic  behaviour  of  a 
thin  resin  layer.  Satisfactory  results  were  obtained,  however 
when  energy  release  rate  calculations  were  performed  it  was  found 
that  the  strain  energy  was  incorrectly  calculated  for  the  springs. 
It  was  therefore  necessary  to  use  the  forces  and  displacements  to 
calculate  the  energy,  and  this  was  incorporated  in  a  post 
processing  routine. 
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6.  Compressive  Failure  Mechanisms 

The  compressive  strength  of  unidirectional  composites  depends  on 
the  shear  stiffness  of  the  composite  which  prevents  buckling  of 
the  fibres.  Another  key  parameter  is  fibre  misalignment.  This 
causes  shear  stresses  and  hence  shear  strains  which  result  In  an 
increase  in  misalignment.  At  a  certain  critical  compressive 
stress  this  process  becomes  unstable  and  failure  occurs.  ABAQUS 
was  used  to  simulate  this  failure  mecnanlsm  [4].  Initially  it  was 
thought  that  orthotropic  material  properties  could  be  used  with 
the  non-linear  geometry  option.  However  this  was  found  not  to 
work  because  the  material  axes  only  rotate  with  rigid  body 
rotation,  but  not  with  shear  stirain.  An  alternative  approach  was 
used  based  on  REBARs  to  model  the  fibre  stiffness.  These  are 
forced  to  maintain  their  position  relative  to  the  element  and  so 
the  fibre  direction  effectively  changes  with  both  rigid  body 
rotation  and  shear  strain,  as  desired. 

The  material  properties  of  the  continuum  elements  containing  the 
REBARS  “ere  chosen  to  represent  the  transverse  and  shear 
properties  of  the  composite.  The  non-linear  shear  response  is 
critical,  and  so  this  was  modelled  using  a  UMAT  subroutine,  as 
before. 

This  modelling  approach  worked  satisfactorily,  and  ABAQUS  was  able 
to  simulate  the  instability  leading  to  compressive  failure.  The 
effect  of  misalignment  angle  on  compressive  strength  was 
investigated  and  shown  to  be  crucial  when  the  composite  is  subject 
to  a  uniform  stress  field  with  no  constraint.  However  when  the 
constraint  imposed  by  the  test  fixture  is  included,  the  effect  of 
misalignment  is  greatly  reduced. 


7.  Micromechanical  Modelling  of  Metal  Matrix  Composites 

A  slice  through  a  section  of  unidirectional  silicon 
carbide-aluminium  was  modelled  using  generalized  plane  strain 
elements  [5].  These  allow  a  constant  but  non-zero  strain  in  the 
fibre  direction.  A  quarter  fibre  with  surrounding  matrix  was 
modelled,  with  boundary  conditions  to  represent  an  infinitely 
repeating  array  of  fibres.  Elastic-plastic  matrix  properties  were 
used,  with  a  Von  Mises  yield  criterion.  Initially  Isotropic 
hardening  was  used,  but  later  the  kinematic  hardening  option  was 
adopted  as  this  is  a  better  model  when  reversed  loading  is 
involved.  The  fibres  were  assumed  to  be  elastic.  Initially  the 
cooldown  from  manufacturing  was  analysed.  This  requires 
temperature  dependent  material  properties  to  be  defined,  which  can 
be  done  in  ABAQUS.  Later  it  was  found  that  this  was  unnecessary, 
and  all  that  was  needed  was  to  apply  the  total  differential 
thermal  strain  and  use  the  room  temperature  properties.  Before 
analysing  the  composite  under  load,  it  was  necessary  to  alter  the 
matrix  material  properties  to  account  for  the  age  hardening  that 
takes  place  after  manufacture.  This  could  be  done  in  ABAQUS 
during  the  analysis  by  using  the  *FIELD  option  or  by  defining  a 
fictitious  temperature  dependence  of  properties. 

A  user  defined  element  was  written  using  the  UEL  subroutine  in 
order  to  model  the  interface.  This  allowed  interface  failure  to 
be  modelled  using  an  interaction  equation  between  normal  and 
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tangential  stresses.  It  also  allowed  a  special  solution  procedure 
to  be  implemented  to  overcome  convergence  problems  during  unstable 
interface  debonding  [6]. 

Transverse  tensile  loading  of  the  composite  was  then  analysed  and 
factors  affecting  the  strength  investigated.  It  was  concluded 
that  the  most  critical  parameters  are  the  Interface  strength  and 
the  level  of  residual  stresses.  The  latter  are  beneficial  to 
transverse  strength  because  they  are  compressive  across  the 
interface,  effectively  increasing  the  interface  strength. 


8.  Concluding  Remarks 

ABAQUS  has  been  used  in  a  wide  range  of  applications  and  has  been 
found  to  be  very  good  for  modelling  composites.  It  is  especially 
suitable  for  non-linear  analysis.  A  good  range  of  standard 
capabilities  is  available,  and  these  can  be  extended  by  means  of 
user  written  subroutines  when  necessary.  This  offers  the  best  of 
both  worlds  in  enabling  the  user  to  benefit  from  a  powerful, 
general  purpose,  commercially  supported  program  whilst  being  able 
to  implement  specific  capabilities  where  required  via  user 
subroutines.  Whilst  ABAQUS  has  limitations,  and  some  problems 
have  been  encountered,  it  is  a  program  undergoing  rapid 
development,  and  the  suppliers,  Hibbitt,  Karlsson  and  Sorensen 
Inc.,  have  generally  been  found  to  be  responsive  to  problems  that 
have  arisen. 
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ANALYSIS  OF  COMPOSITE  STRIJCTURFS 


By  W.S.Amold,  I.H.Marshall,  J.M.Thompson  and  J.Wood 
Deparuneni  of  Mechanical  and  Manufaciuhng  Engineering 
Paisley  College 

1.  INTRODUCnON 

Due  10  certain  intrinsic  fa  tors  associated  with  composite  materials  (  for  example  their  laminar  onhotropic  na¬ 
ture  ),  many  problems  are  extremely  difficult  to  study  using  classical  techniques  or  are  completely  intractable 
with  no  closed  form  solu'ions  available.  Using  Finite  Elements  to  analyse  such  problems  can  often  appear  an 
attractive  alternative,  however,  there  are  several  distinctions  which  have  to  be  acknowledged  when  using  the 
technique  in  comparison  to  the  analysis  of  homogenous  isotropic  matenals. 

This  paper  outlines  some  of  these  considerations  with  reference  to  element  types  and  modelling  strategies  avail¬ 
able  in  commercial  F  E.  packages.  Some  inibal  thoughts  on  benchmarks  for  composite  elements  are  also  pres¬ 
ented. 

2.  ELEMENT  TYPES 

As  with  conventional  matenals  the  correct  choice  of  element  type  is  cnucal.  Clearly,  the  element  formulauon 
should  be  appropnate  to  the  type  of  problem  and  the  analyst  should  be  aware  of  its  limitauons  in  terms  of  dis¬ 
placements  and  in  particular  stress  output.  Four  broad  categories  are  available  in  commercial  systems. 

(1)  Laminated  plaie/shell  elements  ;  These  elements  are  used  for  the  Represcniauve  Volume  Element  tKVE) 
modelling  of  laminated  plates  and  shells. 

(2)  Sandwich  plate/shell  elements  ;  These  elements  are  used  to  model  'sandwich'  type  structures  such  as  foam 
filled  plastics  and  honeycomb  metallic  structures. 

(3)  Orthotropic  brick  :  This  element  could  be  used  for  RVE  modelling  of  a  laminate  or  ply  by  ply  modelling  in 
a  Replica  analysis.  However,  the  element  would  not  allow  individual  ply  stresses  to  be  obtained  in  an  RVE 
model,  consequently  its  main  use  is  in  Replica  modelling,  however  the  element  is  widely  used  for  structural 
stiffness  assessment  and  dynamic  analysis. 

(4)  Laminated  or  stacked  bnck  :  This  element  is  used  for  the  RVE  modelling  of  laminated  composites.  Inter¬ 
laminar  shear  and  normal  stresses  are  normally  available  as  output.  However  the  accuracy  of  these  quanuties 
shall  be  a  function  of  the  number  of  plies  per  element. 


3.  MODELLING  TECHNIQimS 

Various  levels  of  analysis  are  possible,  from  the  microscopic  level  in  which  the  composite  consdtuents  are  ana¬ 
lysed  to  a  ply  by  ply  analysis  in  which  the  plies  are  modelled  as  distinct  enbiies  through  to  a  macro  analysis 
where  equivalent  properties  are  evaluated  for  the  laminate,  it  is  then  treated  as  an  onhotropic  homogeneous  ma¬ 
terial. 

The  type  of  model  required  for  a  particular  analysis  is  very  much  dependant  on  the  engineering  requirements. 
Various  modelling  techniques  uulised  in  the  analysis  of  composite  structures  are  shown  in  Figure  1. 
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(1)  Represeniative  Volume  Element  (RVE)  Models  ;  This  approach  is  based  on  the  assumption  that  the  stress 
response  of  a  composite  can  be  determined  by  the  volume  average  response  of  an  element.  The  properties  of  the 
RVE  element  at  the  macro  level  are  derived  from  individual  fibre/mairix/ply  properties,  using  lamination 
theory,  the  rule  of  mixtures  or  some  other  micro  mechanics  theory.  Alternatively,  these  properties  can  be  deter¬ 
mined  from  various  tests  performed  on  representative  specimens  of  the  material. 

The  RVE  approach  generally  simplifies  the  analysis  problem  and  allows  ply  stresses  to  be  obtained  as  a  post 
processing  operation  from  laminate  strains.  However  this  approach  is  not  generally  employed  to  analyse  de¬ 
tailed  conditions  such  as  'free  edge  effects',  damage  modelling  or  geometrical  discontinuities  such  as  ply  drop 
offs. 

A  numacr  of  'hybnd'  RVE  elements  have  been  derived  which  allows  the  stress  singularity  at  free  edges  to  be 
modelid.  However  none  of  these  appear  to  have  been  implemented  in  commercial  finite  element  systems. 

(2)  Ri  plica  Modelling  (3-D  Modelling)  ;  In  this  approach,  laminate  plies  or  in  some  cases  discrete  fibre/mairix 
models  are  represented  by  a  series  of  unlaminatcd  orthotropic  elements,  perhaps  as  many  as  two  or  three  para¬ 
bolic  elements  per  ply.  Obviously,  this  level  of  analysis  cannot  be  used  for  large  snuctural  problems  due  to 
computational  limitations  but  is  of  value  when  looking  at  detailed  interlaminar  spess  distnbuDons  and  damage 
mechanics. 

(3)  Hybnd  Modelling  ;  With  this  approach,  volume  average  property  distnbution  modelling  (RVE)  is  combuied 
with  discrete  ply  modelling  (Replica)  in  the  same  analysis  to  create  global/tocal  models.  Care  must  be  taken  in 
the  selection  of  the  global/local  interface  assumptions.  With  this  approach,  stacked  bricks  or  laminated  plate 
platc/shell  elements  can  be  used  in  conjunction  with  a  local  ply  by  ply  model. 

(4)  Nested  Modelling  ;  This  is  exactly  the  same  approach  adopted  traditionally  for  the  analysis  of  detail  in  iso- 
uropic  structures.  An  area  of  a  coarsely  meshed  model ’s  isolated  and  re-meshed  with  either  a  finer  mesh  of 
RVE  elements  or  a  new  mesh  of  Replica  elements.  The  boundary  conditions  applied  to  the  nested  model  are 
obtained  from  the  results  corresponding  to  the  boundary  region  in  the  coarse  model. 

The  most  managable  level  of  analysis  involves  the  use  of  plate/shell  elements  to  model  the  'cwo  dimensional' 
response  of  structures.  This  type  of  analysis  can  be  used  to  predict  the  overall  structural  stiffness  and  in-plane 
stress  distribution  but  ui  most  cases  not  its  failure.  This  is  because  failures  tend  to  be  uiibaied  at  discontinuities 
which  require  a  lower  level  of  analysis.  However,  it  is  worth  menuoning  that  some  semi  -  empirical  methods 
have  been  developed  to  predict  failure  at  SCFs  which  use  two-dimensional  '  elements,  e  g.  the  Charactensuc 
Distance  approach  by  Whitney  and  Nuismer  [I). 

4.  OTHER  CONSIDERATIONS 

Oassical  lamination  theory,  which  is  based  on  a  plane  stress  assumpuon,  cannot  lepresent  Lhe  uiierlaminar 
shear  and  normal  stress  components.  Although  these  stresses  may  be  small  in  relation  to  the  ui-plane  stresses, 
the  strengths  associated  with  them  are  also  invariably  much  smaller.  These  components  of  stress  are  therefore 
extremely  important  in  causmg  the  failures  which  are  unique  to  composite  materials. 

Normally,  the  only  way  of  accurately  obtaining  Jiese  quantities  is  to  use  a  ply  by  ply  Replica  model  or  a  hybrid 
or  nested  approach.  However  most  commercial  finite  element  systems  will  attempt  to  approximate  the  inter¬ 
laminar  shear  stress  by  manipulaung  the  resulting  displacements  and  strains. 

For  homogeneous  isotropic  materials  the  applicauon  limits  ot  Kunhoff  plate  theory  (  transverse  shear  effects 
igiKired  )  are  relatively  well  known  and  applicable  to  thin  plates.  For  thicker  plates  the  inclusion  of  Uansverse 
shear  effects  is  included,  e.g.  Mindlin  type  elements  which  assume  constant  transverse  shear  through  the  thick¬ 
ness.  Very  thick  plates  must  be  modeled  using  three  dimensional  bnek  elements. 
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Composite  materials  however,  exhibit  EAj  rauos  which  may  be  an  order  of  magnitude  greater  than  in  the  iso¬ 
tropic  case.  With  a  relauvely  low  shear  modulus,  it  is  known  that  the  level  of  transverse  shear  deformauon  is 
important  in  determining  maximum  deflecuons.  vibration  natural  frerjuencies,  buckling  loads  and  suess  results 
(2].  In  most  cases,  plates  which  would  normally  be  caiegonzed  as  being  'thin'  require  transverse  shearing  ef¬ 
fects  to  be  included  for  accurate  results. 

S.  BENCHMARKS 

With  the  many  claims  made  by  software  vendors  concerning  element  performance  it  is  important,  as  a  matter  of 
quality  assurance  that  valid  benchmarks  for  composites  are  established.  Agencies  such  as  NAFEMS  are  current¬ 
ly  addressing  this  issue  (3).  There  are  several  desirable  features  which  benchmarks  should  possess,  viz  -  (a)  as 
simple  as  possible  (b)  only  one  variable  at  a  time  should  be  considered  (c)  consider  only  clearly  valid  targets 
either  from  theory  or  from  converged  Replica  modelling. 

Furthermore  there  are  several  generic  features  which  benchmarks  should  address,  including  ; 

-  effect  of  lack  of  inter-element  fibre  continuity  as  a  results  of  the  element  material  model  or  as  a  result  of  geo¬ 
metrical  approximation. 

-  membraneAiending/twisting  coupling  for  unsymmetrical  layups 

-  ability  to  model  SCF's  in  orlhotropic  materials  due  to  the  high  stress  gradient  piesent  in  orihotropic  struemres 

-  in  general,  composite  benchmarks  should  conform  to  isotropic  benchmarks. 

-  the  recovery  of  interlaminar  shear  and  normal  stresses  particularly  for  RVE  models. 

-  examine  the  effect  of  material  aspect  ratios 

Clearly,  other  important  areas  should  be  addressed,  including  dynamic,  thermal  and  non  -linear  matcnal  mo¬ 
dels.  Two  simple  examples  are  now  considered  which  highlight  some  of  the  above  features. 

EXAMPLE  A  .VfEMBRANE  RESPONSE  QF  COMPOSITE  CYLINDERS 

As  shown  in  figure  2,  the  problem  consists  of  an  open  ended  cybnder  under  internal  pressure.  Two  models  are 
considered.  In  the  first,  a  quarter  model  is  used  to  represent  an  orthotropic  cylinder  whose  pnncipal  material 
direction  is  in  the  hoop  direction  and  is  coincident  with  the  element  axis.  In  the  second  case,  a  wound  cylinder 
in  which  the  principal  material  direction  is  rotated  from  the  clement  axis  is  represented  by  a  full  model.  The 
analysis  was  carried  out  using  P/FEA  (4).  Two  element  types  were  considered,  the  QUAD/8  doubly  curved 
shell  with  transverse  shear  effects  and  the  linear  QUAD/4  shell. 

In  both  models  only  one  geometrical  patch  was  used  every  90  degree  segment  to  generate  the  elements  thereby 
avoiding  any  disparity  due  to  geometrical  inaccuracies.  Also  u  II  cases  the  element  aspect  rauos  were  held 
constant 

Ca.se  1 

The  object  of  the  first  model  is  to  establish  the  sensiuvity  of  the  hoop  stresses  to  the  number  of  elements  sub¬ 
tending  the  90  degree  section  for  both  an  isotropic  and  orthouopic  material.  Results  are  shown  in  Table  1. 

Only  two  QUAD/8  elements  are  required  to  reach  the  target  value  while  the  QUAD/4  results  converge  towards 
the  target  with  eight  elements.  Both  material  models  give  idenucal  results  with  each  element  type  indicating 
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that  there  is  no  additional  effect  of  inter-element  fibre  continuity  m  the  case  where  the  principal  material  direc¬ 
tion  is  coincident  with  the  element  axis.  The  poorer  results  of  the  linear  elements  were  expected  due  to  the  seg¬ 
mented  representation  of  the  cylindrical  surface.  This  is  designated  method  1. 

Case  2 

The  object  of  this  model  is  to  determine  the  sensitivity  of  the  material  stresses  to  the  number  of  elements  sub¬ 
tending  90  degrees  for  the  case  where  the  principal  material  direction  is  rotated  with  respect  to  the  element  axis. 
This  is  designated  method  2. 

Two  methrxis  are  available  in  the  Patran  to  model  material  directions.  In  the  fust,  the  material  can  be  defined 
as  a  one  ply  laminate  in  which  the  principal  laminate  direction  is  coincident  with  the  element  axis  i.e.  the  hoop 
direction.  The  laminate  is  then  defined  to  be  composed  of  one  ply  orientated  at  an  angle  of  60  degrees  to  the 
element  axis. 

In  the  second  method  the  principal  laminate  direction  is  defined  as  lying  at  60  degrees  to  the  element  axis  .  The 
laminate  is  then  defined  as  being  composed  of  one  ply  orientated  at  0  degrees  to  the  laminate  direction.  Results 
are  shown  in  table  2. 

Using  the  QUAD/8  element  and  method  1  to  generate  the  material  directions,  gives  the  target  values  for  model 
WCYL2LO.  Employing  method  2  however,  gives  errors  of  14%  and  5%  on  the  mauirial  stresses.  Results  im¬ 
prove  significantly  by  doubling  the  number  of  elements  around  the  circumference.  Results  for  the  QUAD/4  el¬ 
ements  again  show  disparities  between  the  two  methods.  Doubling  the  number  of  elements  gives  only  minor 
improvements  in  stress  with  respect  to  the  target  values. 

These  apparent  dispcrepancies  are  currently  being  investigated  by  the  software  suppliers.  This  example  serves 
to  illustrate  that  simple  checks  of  this  nature  can  indicate  basic  innacutacies  in  the  modelling  of  material  prop¬ 
erties.  Bearing  in  mind  the  low  transverse  strengths  of  composites,  even  small  errors  in  stress  can  be  important 
Additionally,  in  this  example  the  cylinder  is  essentially  subjected  to  a  uniaxial  stress  field.  The  previous  dispa¬ 
rities  may  be  more  acute  in  shells  subjected  to  bending  or  at  SCF’s. 

EXAMPLE  B  BENDING  OF  A  TWO  LAYERED  COMPOSITE  BEAM 

This  example  consists  of  a  two  layered  cantilever  beam  subjected  to  an  end  load  and  applied  couple  as  shown  in 
Figures  3  and  4.  The  problem  has  been  solved  by  Lekhinitskii  (5)  who  provides  bending  stress  and  shear  stress 
distributions  through  the  beam  thickness  thereby  enabling  direct  comparisons  with  F.E.  solutions. 

As  the  theoretical  model  is  based  on  Kirkhoff  theory,  the  restraints  do  not  fully  fix  the  end  nodes  but  allow  for 
transverse  deformation  relative  to  the  neutral  axis.  The  ratio  of  Young’s  Moduli  in  the  longitudinal  direction  of 
the  beam  for  the  two  layers  is  9E1  =  E2.  The  analysis  was  carried  out  using  PAFEC. 

Results  are  ploaed  in  Figure  S  for  the  bending  stress  and  interlaminar  shear  stress.  As  can  be  seen  the  results 
are  in  excellent  agreement  The  high  mesh  refinement  used  in  this  model  is  required  for  accurate  retrieval  of  the 
interlaminar  shear  stresses.  Accurate  bending  stresses  were  recovered  using  models  with  only  three  elements 
through  the  thickness. 

It  is  considered  that  this  example  would  make  an  excellent  benchmark  for  Replica  modelling  problems  where 
the  accuracy  of  bending  and  interlaminar  shear  stresses  are  important  parameters  in  the  analysis. 

6.  CONCLUDING  COMMENTS 

Various  modelling  strategies  and  observations  on  the  importance  of  shear  deformation  in  composite  behaviour 
are  indicated.  Important  aspects  on  benchmarks  for  composite  elements  are  discussed  and  two  sample  problems 
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arc  presented. 
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MODEL 

ELEMENT 

EL/90  degree 

MATERIAL 

TARGET 

RESULTS 

CYL2I 

QUAD/8 

2 

ISO 

10 

10 

CYL2AI 

QUAD/4 

2 

ISO 

10 

9.24 

CVTAAI 

QUAD/4 

4 

ISO 

10 

9.81 

CYL6AI 

QUAD/4 

6 

ISO 

10 

9.92 

CYL8AI 

QUAD/4 

8 

ISO 

10 

9.96 

CYL20 

QUAD/8 

2 

ORTH 

10 

10 

CYL2AO 

QUAD/4 

2 

ORTH 

10 

9.24 

CYL8AO 

QUAD/4 

8 

ORTH 

10 

9.96 

QUARTER  MODEL  RESULTS  -  EXAMPLE  A 


TABLE  1 


MODEL 

ELEMENT 

EL/90  degree 

METHOD 

TARGET 

L  T 

RESULTS 

L  T 

WC\X2LO 

QUAD/8 

2 

1 

2.5 

7,5 

2.5 

7.5 

WCYL20 

QUAD/8 

2 

2 

2J 

7.5 

2.14 

7.86 

WC'yL40 

QUAD/8 

4 

2 

2.5 

7.5 

2.41 

7.59 

WCYL4ALO 

QUAD/4 

4 

1 

2.5 

7.5 

2.38 

7.4 

WCYL4AO 

QUAD/4 

4 

2 

2.5 

7.5 

2.48 

7.32 

WCYL8AO 

QUAD/4 

8 

2 

2.5 

7.5 

2.41 

7.42 

FULL  MODEL  RESULTS-  EXAMPLE  A 


TABLE  2 
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FINITE  ELEMENT  MODELLING 
OF  COMPOSITE  MATERIALS  AND  STRUCTURES 


REPRESENTATIVE  VOLUME  REPLICA  MODELLING 

ELEMENTS  (  RVE’s ) 


-  rule  of  mixtures 

-  laminadon  theory 

-  shell  elements 

-  laminated  bricks 

-  allows  ply  stresses  to  be 
obtained  as  a  post  -processing 
operation 

-  validity  ? 

-  accuracy  ? 


-  orthocropic  bricks 

-  ply  by  ply  models 

-  3-4  elements  per  ply 

-  post  processing  cumbersome 

-  extremely  large  models 

-  local  and  global  effects  theoetically  possible 


HYBRID  MODELLING  or  NESTED  MODELLING 


RVE  followed  by  local 
REPLICA  modelling 

-  combinations  of  laminated  shell, 
bricks  and  orthotropib  bricks 


RVE  and  REPLICA 
in  same  model 


HGURE  1 
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CASE  I 


PRESSURE  T 

X  ^  (3) 


[  1  rcsiiaints 

CYLINDER  DIMENSIONS  ; 
I  =  1  mm 
r  =  10  mm 
L  =  100  mm 

Internal  pressure  =  1  N/mm2 


QUARTER  MODEL 


HELIX  ANGLE  =  30  degrees 


EXAMPLE  A  MEMBRANE  RESPONSE  OF  COMPOSITE  CYLINDERS  RGURE  2 
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Finite  element  O 


xtOO 


(at  itz3QI 


tat  t  s  30) 


EXAMPLE  B  BENDING  STRESS  AND  SHEAR  STRESS  DISTRIBUTION 

RGURE  5 
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Finite  Element  Analysis  of  Composite  Space  Systems 


Ml  Gunn  and  LHollawav 
Composite  structures  tor  space 

For  a  number  of  years  the  Composite  Structures  Research  Unit  in  the 
Department  of  Civil  fngineering  at  the  University  of  Surrey  has  been 
concerned  with  the  design  of  large  structures  for  use  in  space.  This  type  of 
structure  may  be  used  for  space  communications,  space  science  or  the 
establishment  of  the  first  permanent  space  station.  Composite  materials  are 
particularly  appropriate  for  this  typo  of  application  because  of  their 
favourable  strength  and  stiffness  to  mass  ratios.  A  typical  structure  is  a  three 
dimensional  double  layer  truss  designed  to  support  a  parabolic  reflector  for 
communications.  These  structures  present  a  number  of  design  problems 
where  finite  element  analysis  may  be  of  assistance,  including  vibrational 
characteristics,  deployment  from  their  initial  (folded)  form,  and  the 
effectiveness  of  the  jointing  system  that  is  adopted. 

Using  more  than  one  program 

The  analyses  reported  below  have  taken  place  over  the  last  six  years. 
Initially  we  only  had  access  to  LUSAS  at  the  University  of  Surrey.  At  the 
time  of  writing  we  have  LUSAS,  ABACUS,  NASTRAN  and  ANSYS  available 
for  use.  There  are  a  number  of  good  reasons  for  using  more  than  one  finite 
element  program  in  research  projects.  Firstly  programs  vary  in  their 
capabilities  and  ease  of  use.  Secondly,  there  is  the  question  of  reliability  of 
results.  Repeating  an  analysis  using  a  different  program  checks  data 
perpetration  and  understanding  of  how  to  operate  program  options  as  well 
as  the  correctness  of  the  program  itself. 

The  first  two  analyses  would  seem  to  demonstrate  a  general  superiority  of 
ABACUS  over  LUSAS.  We  should  empasise  that  we  have  faith  in  the 
LUSAS  program  and  continue  to  use  it  in  our  research.  A  significant 
advantage  of  LUSAS  is  that  it  is  a  British  program  and  excellent  local 
support  is  available.  ABACUS  is  an  American  program  and  although  it  is 
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probably  the  leading  non  linear  code  available  one  can  expect  difficulties  in 
obtaining  the  same  level  of  support  as  for  a  program  originating  in  the  UK. 

Dynamic  analysis  of  unsupported  structures 

It  is  well  known  that  a  structure  must  have  a  minimum  number  of  supports  if 
one  is  to  perform  a  proper  static  analysis.  A  structure  in  space  is 
unsupported  and  therefore  a  step  by  step  dynamic  analysis  (based  on 
Newton's  laws  of  motion  rather  than  static  equilibrium)  is  appropriate.  The 
absence  of  the  section  of  data  that  describes  the  support  conditions  led  to  a 
data  error  from  LUSAS.  The  insertion  of  "free"  boundary  conditions  led  to  a 
division  by  zero  error.  The  question  then  arose;  was  this  a  problem  with 
LUSAS  or  Newmark's  method?  The  various  step  by  step  dynamic 
algorithms  described  by  Bathe  (1977)  were  programmed  and  tested  against 
a  simple  problem  with  an  analytical  solution.  None  had  a  problem  with  the 
absence  of  displacement  controlled  boundary  conditions.  In  the  meantime 
ABACUS  was  discovered  mounted  at  the  University  of  Manchester  Regional 
Computer  Centre.  One  of  the  problems  in  the  Examples  Manual  was  the 
step  by  step  dynamic  analysis  of  an  unrestrained  stmcture.  ABACUS  was 
then  successfully  applied  to  our  problem  (Hollaway  et  al,  1987).  A 
particularly  useful  feature  of  ABACUS  was  the  possibility  of  automatically 
determining  the  size  of  time  step  required  as  the  analysis  progressed.  This 
is  done  by  estimating  the  equilibrium  error  at  the  mid  point  of  the  time  step 
and  repeating  the  step  (with  a  shortened  time)  if  this  error  exceeds  some 
specified  value. 

Stress  analysis  of  a  joint 

The  second  problem  to  be  discussed  is  the  non  linear  analysis  of  a  tubular 
lap  joint.  Various  analyses  were  carried  out  with  tubes  made  from  either 
aluminium  or  carbon  fibre/polyethersulphone  (PES)  composite  (Romhi  el  al, 
1986,  Romhi,  1990).  One  variation  on  the  basic  configuration  was  that  the 
outer  tube  was  replaced  by  an  end  cap  which  would  enable  the  member  to 
be  joined  to  others  in  a  three  dimensional  truss  structure.  The  objective  of 
the  analyses  was  to  examine  the  stress  distribution  in  the  joint  and  assess 
the  effect  of  the  different  material  properties  on  the  behaviour  of  the  joint  as 
failure  was  approached.  Axial  loads  were  applied  to  the  joint  and  the 
condition  analysed  was  one  of  axisymmetric  plane  strain.  The  tubes  were 
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assigned  elastic  properties.  The  adhesive  was  modelled  as  an 
elastoplasitic  material,  yielding  according  to  Von  Mises  yield  condition. 

Both  LUSAS  and  ABACUS  showed  yielding  starting  at  the  ends  of  the 
adhesive.  Under  increased  load  these  zones  of  yielding  spread  towards  the 
centre  of  the  joint.  At  the  time  these  analyses  were  performed  ABACUS  had 
an  automatic  load  incrementation  facility  similar  to  that  described  above  for 
dynamic  analysis  whereas  LUSAS  required  that  the  size  of  each  load  step 
be  specified  before  the  analysis  started.  As  a  result  we  had  difficulties  in 
obtaining  a  converged  analysis  using  LUSAS  despite  several  repetions  with 
modified  sizes  of  load  increments. 

Using  the  "User  Material"  facility  in  ABACUS  we  were  able  to  modify  the 
yield  condition  for  the  adhesive  to  one  where  yielding  depended  on  th« 
mean  normal  pressure.  This  involved  producing  a  Fortran  subroutine  which 
would  calculate  the  elastoplastic  modulus  matrix  relating  an  increment  of 
strain  to  an  increment  of  stress  and  in  addition  would  ensure  that  stresses 
did  not  exceed  yield  by  performing  a  plastic  stress  correction.  We  should 
issue  a  warning  that  using  this  type  of  facility  requires  expert  knowledge  of 
finite  element  programming.  The  work  was  checked  by  programming 
ABACUS  with  Von  Mises  yield  condition  via  the  User  Material  facility  and 
checking  against  the  standard  Von  Mises  code  in  ABACUS.  All  the  routines 
that  we  wrote  were  also  independently  checked  in  another  (less 
sophisticated)  finite  element  program. 

Specification  of  anisotropic  properties 

Another  research  project  examined  the  behaviour  of  a  single  composite  tube 
under  impact  and  general  vibrational  loading.  Here  the  analysis  was  to  be 
three  dimensional  and  the  first  question  considered  was  whether  the  tube 
should  be  modelled  by  this  or  thick  shell  elements.  The  research  worker 
had  taken  care  to  establish  the  anisotropic  elastic  properties  of  the  carbon 
and  glass  reinforced  PES  composites,  but  there  were  unexpected 
difficulties  in  using  these  properties  in  the  finite  element  models.  The  most 
straight  forward  case  was  where  the  semiloof  curved  thin  shell  elements  in 
LUSAS  were  used.  Properties  for  these  elements  are  always  specified  in 
terms  of  local  co-ordinates  associated  with  each  element.  The  thick  shell 
elements  in  LUSAS  are  basically  "cut-down*  versions  of  the  familiar  20 
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noded  isoparametric  brick  elements  with  the  edge  nodes  removed  through 
the  thickness  of  shell.  By  default  properties  are  assumed  to  be  in  terms  of 
the  global  axis  system.  To  enter  properties  aligned  with  the  element 
geometry  it  is  necessary  to  specify  a  local  cartesian  coordinate  system  for 
each  element.  This  is  unsatisfactory  for  two  reasons.  Firstly  an  extra 
approximation  is  introduced  into  the  analysis  as  the  curvature  of  the  element 
is  ignored.  Secondly  it  is  strictly  speaking,  unnecessary  since  this 
information  has  already  been  entered  in  the  element  connectivity 
information. 

The  use  of  ABAQUS  was  considered,  but  not  pursued,  as  properties  had  to 
be  entered  directly  as  elements  of  the  tensor  Dp.  Finally  a  reduced  from  of 
anisotropy  was  used  with  LUSAS  in  conjunction  with  properties  specified  in 
the  directions  of  the  global  axes.  One  global  axis  was  aligned  with  the  axis 
of  the  tube  (the  90%  fibre  direction)  and  all  properties  in  perpendicular 
planes  were  assumed  to  be  equal.  In  practice,  the  results  of  the  analyses 
conducted  were  mainly  controlled  by  the  longitudinal  modulus  and 
satisfactory  modelling  was  achieved  with  isotropic  properties.  Hence  in  this 
case  the  agonizing  over  correctly  entering  anisotropic  properties  was 
unnecessary. 

Analysis  of  a  plate  loaded  bv  a  single  pin 

This  is  the  simplest  form  of  analysis  relevant  to  bolted  joint  behaviour.  This 
problem  was  used  to  compare  ease  of  data  generation,  modelling  options 
and  accuracy  of  calculated  stresses  for  LUSAS,  ABACUS  and 
MSC/NASTFIAN.  All  three  packages  have  internal  geometric  preprocessors 
and  in  summary:  LUSAS's  preprocessor  is  simple  to  use  but  less  general 
and  NASTRAN's  preprocessor  has  the  most  generality  but  suffers  from 
complexity. 

ABAQUS  and  NASTRAN  both  have  options  for  calculating  the  effective 
material  properties  according  to  laminated  plate  theory  and  can  also  output 
stresses  in  the  different  layers.  LUSAS  and  ABAQUS  have  similar  options 
for  printing  stresses  at  nodal  points:  either  averaged  or  based  on  individual 
element  values.  MSC/NASTRAN,  on  the  other  hand,  calculates  stresses  at 
the  element  corners  and  centre  expressed  in  the  element's  local  co¬ 
ordinate  system.  The  original  intention  had  been  to  use  the  8-noded 
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isoparametric  element  with  reduced  integration  in  each  package. 
Unfortunately,  NASTRAN  does  not  allow  a  global  definition  of  the  material 
axes  with  this  element  type  and  so  a  finer  mesh  of  4-noded  elements  was 
used  instead.  Analyses  were  also  conducted  with  9-noded  elements  and 
the  8-noded  elements  with  full  integration. 

The  results  were  assessed  in  terms  of  the  overall  stiffness  of  the  plate  (all 
element  types  and  packages  gave  similar  answers  and  were  considered 
adequate)  and  the  stress  distribution  around  the  hole.  Overall  the  best 
results  were  obtained  using  the  9-noded  elements  and  the  8-noded 
elements  with  full  integration,  although  all  element  types  gave  similar 
answers  for  hoop  and  radial  stresses  around  the  pin.  There  were  significant 
differences  in  the  shear  stress  distributions  calculated  with  the  worst  results 
obtained  form  the  4-noded  elements. 
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FINITE  ELEMENT  ANALYSIS  OF  COMPOSITE  NAVAL  STRUCTURES 


J  D  McVEE  and  R  S  DOW 
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Dunfermline,  Fife,  KYll  5PW 

An  overview  is  presented  of  the  application  of  finite 
element  techniques  to  the  structural  analysis  and  design  of 
composite  naval  structures  carried  out  during  the  last  twenty- 
five  years  (1)*  at  the  Dunfermline  site  of  the  recently  formed 
Defence  Research  Agency,  Maritime  Division.  It  includes  a  brief 
description  of  marine  applications  of  laminated  glass  reinforced 
plastic  (GRP)  composites  along  with  a  short  discussion  of  typical 
laminate  properties,  permissible  stress  levels  and  relevant 
failure  mechanisms.  Where  applicable,  other  theoretical  methods 
such  as  laminated  orthotropic  beam  theory,  folded  plate 
techniques  and  axisymmetric  shell  stress  and  buckling  analysis 
are  also  considered.  The  importance  of  theoretical  results 
correlation  with  full  and  model  scale  test  data  is  emphasised. 

In  particular,  extensive  use  of  several  linear  and  non¬ 
linear  proprietary  finite  element  computer  codes  in  the  evolution 
(2,  3,  4)  of  three  classes  of  GRP  mine  countermeasures  vessels 
(MCMVs)  developed  by  the  Royal  Navy  during  the  last  thirty  years 
is  very  briefly  described.  Examples  are  given  of  primarily 
linear  finite  element  analysis  ranging  from  the  determination  of 
the  global  elastic  response  of  large  scale  MCMV  test  sections, 
corrugated  GRP  hulls,  and  subsequent  detailed  examination  of 
various  GRP  joint  designs  to  the  prediction  of  stress 
distributions  in  typical  GRP  interlaminar  tensile  strength 
specimens.  Non-linear  analysis  concerned  with  compressive 
buckling  of  longitudinally  stiffened  panels  (5),  interactive 
buckling  effects  (6)  and  the  initial  buckling  and  post  buckling 
of  grillages  subjected  to  combined  shear  and  direct  loading  (7) 
are  discussed. 

Recent  studies  concerned  with  the  potential  application  of 
hybrid  GRP/steel  superstructures  are  outlined.  These  consist  of 
preliminary  static  elastic  analysis  (8)  which  indicate  that 
fatigue  failure  associated  with  hull-superstructure  interaction 
can  be  significantly  reduced.  Moreover,  dynamic  response 
analysis  (9)  of  representative  full-scale  superstructure  deck 
panels,  and  a  corresponding  three  dimensional  deckhouse  module, 
have  been  correlated  with  experimental  test  data  which 
demonstrates  the  capability  of  this  form  of  construction  to 
withstand  external  air  blast  loading. 

With  reference  to  submarines,  a  series  of  linear  static 
elastic  response  and  buckling  analysis  (10)  performed  on  a 
sequence  of  models  representing  a  full  scale  GRP  deck  casing  test 
section  subjected  to  proof  load  pressures  are  described  which 
confirmed  that  failure  was  due  to  through  thickness  shear 
effects.  Regarding  the  potential  use  of  fibre  reinforced  polymer 
(FRP)  components  in  the  pressure  hulls  of  submersibles  (11),  the 
theoretical  performance  of  FRP  is  contrasted  with  that  of  steel 
and  light  alloy  construction  in  two  illustrative  trial  designs 
corresponding  to  a  shallow  water  manned  submersible  and  a  deep 
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water  autonomous  underwater  vehicle.  Axisymmetric  thin  and  thick 
shell  and  three  dimensional  non-linear  finite  element  discretis¬ 
ations  necessary  to  obtain  critical  buckling  pressures  are 
highl ighted . 

Finally,  some  suggestions  are  made  for  code  improvements 
including  increased  ease  of  data  input,  improved  elements,  more 
general  stress-strain  definitions,  access  to  more  general  failure 
criteria  and  more  user  orientated  visual  pre  and  especially  post 
processing  capabilities. 
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Tel:  071-465-2264 


This  presentation  describes  the  evolution  of  a  method  for  the  analy.'.is  of 
composite  racing  yachts,  along  with  the  development  of  the  required 
software.  It  is  based  on  Arup  experience  in  the  design  and  analysis  of  the 
hull  structure  for  the  International  Offshore  Racing  yacht  'Satquote 
British  Defender'  and  an  America's  Cup  Class  yacht  for  an  Australian 
syndicate.  Most  recently  the  techniques  have  been  used  to  analyse  Dennis 
Conner's  latest  America's  Cup  yacht  'Stars  and  Stripes'. 

The  presentation  discusses: 

'  The  significance  of  the  hull  shell  shape  to  its  structural  behaviour. 

-  How  a  finite  element  model  can  be  created  from  the  available 
geometric  data. 

-  How  material  input  data  can  be  manipulated. 

What  load  cases  need  to  be  considered  and  how  the  loading  is  applied 
to  the  model . 

-  The  special  post  processing  requirements  posed  by  composites  and  what 
data  reductions  techniques  can  be  used. 
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Imperial  College,  London 


This  paper  presents  attempts  and  current  developments  in  softwares  for 
finite  element  aided  design  of  composite  structures.  The  principal  objectives  of 
these  studies  are  to  obtain  methods  and  programs  easy  to  use  and  capable  to  run 
on  micro-computers. 

The  complexity  of  the  behaviour  of  composite  materials  has  led  to  more 
and  more  complex  theories  implemented  in  more  and  more  complex  special 
finite  element  methods.  Presently,  very  sophisticated  and  precise  capabilities  are 
available  for  analysis.  However,  these  softwares  are  less  and  less  easy  to  use  :  the 
collection  and  the  imput  of  the  data  (specially  the  materials  data)  are  long  or 
difficult,  the  selection  of  the  adequate  theory  or  the  adequate  elements  (among  so 
many  available)  may  be  a  real  challenge,  the  programs  generally  require  large 
computers,  finally  the  large  number  of  numerical  results  do  not  receive  easily  a 
precise  and  useful  physicad  meaning. 

On  the  other  hand,  compared  to  these  numerous,  well-developed  and 
complex  methods  for  analysis,  the  design  methods  are  poor.  They  mostly  reduce 
to  very  specific  methods  (generally  using  a  lot  of  algebra)  or  to  optimization 
process  :  using  some  analysis  method  in  each  loop  of  their  iterative  course,  these 
optimization  programs  may  require  huge  computer  facilities  and  long 
computation  times. 

From  these  considerations,  it  appears  useful  to  develop  simple  and  fast 
methods  of  computation.  We  retained  the  features  described  in  the  following. 

1  -  Re-analysis  aided  design  :  Fast  analysis  and  re-analyses  with  finite  elements 
constitute  a  possible  way  for  improving  design  :  instead  of  an  automatic  design 
program,  the  user  has  to  analyse  first,  assess  the  results,  decide  for  the  structural 
changes,  re-analyse  the  structure,  assess  the  new  results,  and  so  on.  A  very  fast 
finite  element  program  is  required  for  that.  We  developed  various 
improvements  in  the  finite  element  method  to  speed  the  computation.  First, 
dimensional  analysis  and  separation  of  the  material  and  the  geometric 
parameters  have  shown  that  pre-computed  quantities  can  be  used  to  form  the 
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stiffness  matrices,  which  results  in  a  faster  formation  of  these  matrices. 
Systematic  use  of  multi-loads,  in  the  initial  analysis,  considerably  shortens  the 
subsequent  re-analyses  for  the  cases  obtained  by  combination  of  solutions.  The 
flexibility  method  proved  to  be  an  efficient  way  for  fast  re-analysis  with  different 
boundary  conditions.  Finally,  perturbation  methods  are  presently  investigated  to 
take  into  account  structural  modifications  in  the  flexibility  matrices. 

2  -  Two-step  search  :  Design  of  composite  laminated  or  sandwich  structures 
should  include  the  choice  of  the  nature  of  the  laminas,  and  the  determination  of 
the  stacking  sequences.  Consequently,  the  complete  problem  of  optimization  -or, 
more  simply,  of  improvement-  of  design  is  ntaterial  and  structural.  To  simplify 
the  approach,  we  introduce  a  two-step  process.  The  first  step  has  to  find  improved 
overall  structural  properties,  while  the  second  step  cortsists  in  finding  a  material 
which  can  reproduce  these  properties.  Two  principal  comments  should  be 
emphasized.  First,  this  two-step  process,  while  simpler,  may  generally  provide  a 
"less  optimal"  solution  than  a  complete  optimization  (if  it  exists  and  can  be 
effectively  processed).  Second,  the  material  design  for  given  structural  properties 
is  questionable  :  it  is  rather  doutfull  that  one  can  obtain  the  wanted  structural 
properties  with  the  available  basic  laminas.  However,  we  solved  this  difficulty  by 
extending  the  range  of  search  possibilities  :  it  is  possible  to  define,  in  a 
quantitative  way,  norms  and  deviations  for  anisotropic  structural  properties,  so 
that,  instead  of  solely  searching  for  materials  which  give  exactly  the  required 
structural  properties,  we  can  search  for  materials  which  give  approximately  these 
properties,  within  a  definite  relative  deviation. 

3  -  User  friendly  programs  :  In  addition  to  fast  processing,  already  discussed,  at 
least  two  requirements  must  be  fulfilled  in  order  to  have  user-friendly  softwares. 
They  must  be  interactive,  specially  with  graphic  facilities,  at  the  input  and  output 
phases.  They  must  run  on  largely  available  and  easy-to-use  computers.  To  meet 
these  requirements,  we  retained  micro-computers  (IBM  PC-AT  or  compatible, 
and  Macintosh  H)  and  low  cost  work-stations.  We  also  included  in  our  programs 
(or  plan  to  include  in  the  programs  currently  in  development)  user-friendly 
facilities,  such  as  use  of  mouse,  pop-up  menus,  graphical  windows,  etc..,  in  the 
way  made  popular  by  the  Apple  softwares. 

Two  softwares  have  been  developed  on  the  above  ideas  and  methods.  The 
program  FEAD-LASP  (Finite  Element  Aided  Design  for  Laminated  And 
Sandwich  Plates)  was  developed  from  parts  of  a  rather  classical  finite  element 
program  for  micro-computer.  It  is  limited  to  rectangular  geometries  (composite 
and  sandwich  panels  in  bending  and  in-plane  deformation)  and  pre-defined 
regular  meshes,  but  takes  into  account  completely  general  boundary  and  load 
conditions.  Generality  is  also  available  for  the  material  properties,  which  include 
bending,  in-plane,  bending-membrane  coupling  and  transverse  shear  stiffnesses 
(extended  laminated  plate  theory,  i.e.  classical  laminated  plate  theory  plus  shear 
effects).  A  16-node  thick  plate  element  with  transverse  shear,  in  which  the 
materials  properties  of  the  laminate  or  sandwich  are  input  using  an  equivalent 
material,  was  retained.  A  faster  method  for  the  formation  and  assembly  of  the 
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stiffness  matrices  of  the  elements  was  derived  and  tested.  It  reduces  the 
corresponding  time  very  significantly,  using  precomputation  and  storage  of 
dimensionless  parts  of  the  stiffness  matrix  elements.  Several  micro-computers 
and  work-stations  were  used  to  test  the  numerical  part  :  it  was  found  that  the 
most  powerful  personal  computers  (such  as  Apple  Macintosh  II  and  IBM  PC-AT) 
as  well  as  the  low  cost  work-stations  (such  as  Apxjllo  DN-3000)  are  well  suited  for 
a  fast  numerical  treatment,  i.e.  oiUy  a  few  minutes  for  the  first  analysis,  and  re¬ 
analyses  in  seconds.  A  user-friendly  input  part  for  the  program  was  developed 
and  makes  possible  the  use  of  the  program  by  designers  with  little  or  no 
knowledge  of  the  internal  structuie  and  methods  of  the  program. 

The  program  FEAD-FLEX  (Finite  Element  Aided  Design  by  Flexibility 
method)  is  currently  in  developr.ient  and  test.  It  already  includes  in  its  features 
the  fast  re-analysis  with  varying  the  boundary  conditions.  The  principle  of  the 
method  is  the  following  :  taking  advantage  of  the  fact  that  the  rigid  body  motions 
are  the  singularities  of  the  equilibrium  system  of  discrete  elastic  structures,  a 
general  solution  can  be  expressed  prior  to  any  assignment  of  displacement 
boundary  conditions,  then  the  solution  for  specified  load  and  boundary 
conditions  can  be  determined  by  solving  an  associated  small-size  linear  system. 
With  this  method,  the  re-analyses  under  various  boundary  conditions  are  made 
easier.  Further,  it  can  simplify  the  solutions  for  contact  problems  and  elastic  crack 
propagations  problems,  which  otherwise  need  long  iteration  methods,  and  be 
useful  for  substructuring  of  large  or  repetitive  structures.  Finally,  perturbation 
methods  for  the  flexibility  matrix  are  currently  implemented  in  this  software. 
They  allow  to  take  into  account  local  modifications  of  the  stiffness  (due  to 
geometric  or  material  changes)  without  coming  back  to  the  beginning  of  the 
process,  and  consequently  reduce  the  computation  time.  We  plan  to  use  this 
method  for  various  problems  in  the  mechanics  of  composite  materials,  such  as 
damage  propagation  analysis,  and  form  optimization. 
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ABSTRACT 

Among  the  applications  we  arc  analyzing  currently,  the  truck  frame  is  the  most  interesting  from 
the  point  of  view  of  problems  and  successes  of  the  application  of  finite  element  method  to  a 
composite  strucnire. 

The  frame,  which  is  the  object  of  the  present  study,  is  composed  of  three  substructures 
(Figure  1). 
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Figure  2.  Dimensions  of  the  truck  frame  in  mm. 

There  arc  two  problems : 

•  Stress  analysis  in  tlie  comers 

*  Stress  analysis  in  the  variable  thickness  areas 

The  variable  thickness  problem  was  solved  by  means  of  8  node  and  linear  displacement  brick. 
But  the  comer  was  an  interesting  area  to  analyze. 

Analysis 

Since  the  comer  is  one  of  the  most  critical  substructutres,  the  following  testings  have  been 
carried  out ; 


Figure  3.  Scheme  of  the  tests  1, 2  and  3. 
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The  material  used  is  ; 

Flangei  :  a  unidirectional  hybrid  carbon  fiber/epoxy  resin  (10%)  and  fiberglass/epoxi  resin 
(90%)  composite  material. 

Web  :  a  unidirectional  fiberglass/epoxy  resin  composite  material. 

In  both  cases,  volumen  fraction  is  0.6. 

The  following  ABAQUS  fifteen  finite  elements  have  been  used  for  the  analysis  or  these  three 
testings : 

S4R  4  node,  reduced  integration,  doubly  curved  shell  with  hourglass  control. 

S8R  8  node,  reduced  integration,  doubly  curved  shell. 

S4R5  4  node,  reduced  integration,  doubly  curved  shell  with  hourglass  control,  using  five 

degrees  of  freedom  per  node. 

S8R5  8  node,  reduced  integration,  doubly  curved  shell  ,using  five  degrees  of  freedom  per 
node. 

S9R5  9  node,  reduced  integration,  doubly  curved  sheU,  using  five  degrees  of  freedom  per 
node. 

C3D8  8  node,  linear  displacement  brick. 

C3D8R  8  node,  linear  displacement,  reduced  integration  brick  including  hourglass  control. 
C3D8H  8  node,  linear  displacement,  constant  pressure  brick. 

C3D8RH  8  node,  linear  displacement,  constant  pressure,  reduced  integration  brick  including 
hourglass  control. 

C3D20  20  node,  quadratic  displacement  brick. 

C3D20R  20  node,  quadratic  displacement,  reduced  integration  brick. 

C3D20H  20  node,  quadratic  displacement,  linear  pressure  brick. 

C3D20RH  20  node,  quadratic  displacement,  linear  pressure,  reduced  integration  brick. 

C3D27  27  node,  cubic  displacement. 

C3D27R  27  node,  cubic  displacement,  reduced  integration  brick. 
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Results 


Test  1.  C3D27  element  is  the  most  accurate  (4  %)  though  results  from  20  and  27  node  element 
are  vey  close  to  the  experimental  data.  The  rest  of  the  elements  give  very  poor  results. 

Test  2.  C3D27  and  C3D27R  elements  give  good  results  (  3  and  3.7  %  error,  respectively )  and 
20  node  elements  present  errors  around  10  %.  C3D8(H)  element  error  is  39  %  and  the  rest  of 
8, 4  and  9  node  elements  present  very  poor  results. 

Test  3.  Results  from  20  and  27  node  element  arc  vey  close  to  the  experimental  data  (  errors 
from  6  to  8  %).  The  rest  of  the  elements  give  very  poor  results  except  C3D8(H),  S4R5  and 
S8R5  (  15  to  20  %  error). 
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NUMERICAL  PROBLEMS  ENCOUNTERED  IN  NONLINEAR  FINITE 
ELEMENT  ANALYSIS  OF  FREE  EDGE  DELAMINATION 


J.C.J.  Schellekens 

Delft  University  of  Technology,  Department  of  Civil  Engineering  / 

TNO  Building  and  Construction  Research 
P.O.  Box  5048, 2600  GA  Delft,  Netherlands 

ABSTRACT 

A  failure  mode  often  encountered  in  structural  composite  laminates  is  delamination,  a 
phenomenon  which  has  great  influence  on  the  structural  integrity.  In  the  last  two  decades 
much  research  has  concentrated  on  the  complex  mechanism  of  delamination.  Due  to  the 
anisotropic  material  properties  and  the  varying  fibre  orientations  each  ply  in  the  laminate 
behaves  differently.  Large  edge  stresses  are  then  necessary  to  preserve  compatibility  of  the 
deformations.  These  large  transverse  normal  and  shear  stresses  are  primarily  responsible 
for  the  initiation  of  delamination.  In  the  literature  various  analytical  models  have  been 
developed  for  the  treatment  of  free  edge  stress  distributions,  e.g.  by  Pagano.  Furthermore 
procedures  to  predict  delamination  onset  and  growth  based  on  the  principle  of  virtual  crack 
extension  have  been  applied  by  Crossman,  Wang  and  O’Brien  amongst  others.  Because  of 
the  stress  singularities  that  exist  at  the  ply  interfaces  near  the  free  edges  and  the  resulting 
mesh-dependence  it  is  commonly  believed  that  the  use  of  stress-based  failure  criteria  does 
not  produce  relevant  results  as  far  as  free  edge  dclamination  is  concerned.  However,  Kim 
and  Soni  indicated  that  an  average  stress  approach  combined  with  an  anisotropic  failure 
criterion  results  in  an  accurate  prediction  of  the  onset  of  delamination.  The  introduction  of 
the  ply  thickness  in  the  determination  of  stresses  is  essential  in  their  approach. 

In  this  contribution  a  nonlinear  finite  element  procedure  is  presented  for  the  prediction  of 
delamination  onset  and  growth.  The  procedure  accounts  for  thermal  effects  since  we  have 
to  take  into  account  the  initial  stresses  in  the  laminate  due  to  the  forming  process.  Although 
the  failure  criterion  is  based  on  stresses  it  is  shown  that,  when  combined  with  a  softening 
type  of  post-crack  response,  a  stress-based  failure  criterion  results  in  a  mesh-objective 
calculation.  The  performance  of  the  method  is  demonstrated  by  means  of  the  analyses  of 
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free  edge  delamination  in  various  graphite-epoxy  laminates  and  the  numerical  results  are 
compared  with  experimenis. 

In  the  examples  twelve-noded  cubic  generalised  plane  strain  elements  with  three 
translational  degrees-of-freedom  in  each  node  have  been  used.  These  elements,  which  are 
assumed  to  remain  elastic  during  the  loading  process,  give  an  accurate  representation  of  the 
stress  concentrations  near  the  free  edges  without  the  need  for  extreme  mesh  refinement.  In 
contrast  to  other  continuum  elements  generalised  plane  strain  elements  are  not  loaded  by 
nodal  forces  or  prescribed  displacements  but  by  prescribed  strains  normal  to  the  element 
plane.  The  individual  plies  are  connected  by  cubic  3D  line  interface  elements  which  are 
well  suited  for  modelling  the  geometric  discontinuity  that  arises  during  delamination, 
which  can  either  be  gradual  (softening  type  behaviour)  or  perfectly  brittle 

In  the  lecture  the  emphasis  will  be  put  on  the  effects  of  mesh  refinement  and  laminate 
thickness  on  the  ultimate  load  capacity  of  the  laminates.  It  will  be  demonstrated  that  the 
inclusion  of  interface  elements  in  the  finite  element  model  results  in  a  proper  description  of 
size-effects  and  mesh-objectivity  due  to  the  softening  type  of  response  after  cracking. 

Since  the  application  of  a  conventional  Gauss  integration  scheme  for  the  assembly  of  the 
stiffness  matrix  of  the  interface  elements  may  result  in  a  poor  element  performance,  the 
effect  of  other  integration  schemes  on  the  clement  behaviour  is  investigated. 

A  major  drawback  of  strain  or  load-controlled  calculations  is  the  fact  that  no  limit  points  in 
the  load  deflection  curve  can  be  passed.  Riks  developed  an  ’arc-length’  method  to 
overcome  these  limitations.  Herein  the  incremental  load-factor  is  constrained  by  the  norm 
of  the  incremental  displacement  vector.  Although  the  arc-length  control  method  has 
proved  to  be  rather  successful,  it  has  been  reported  to  fail  in  situations  of  highly  localised 
failure.  In  these  cases  the  displacement  norm  should  be  determined  considering  only  the 
dominant  degrees  of  freedom.  Therefore  a  modified  version  of  the  standard  arc-length 
method  will  be  discussed  which  results  in  a  stable  solution  process  upon  passing  limit 
points  in  the  structural  response. 
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Thermal  Analysis  of  Laminated  Composites 
-  Element  Development  and  Experience  with  MSC/NASTRAN 
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Workshop:  Finite  Element  Methods  for  Composites.  September  Sth-6th,  1991,  London 


1  INTRODUCTION 

Thin-walled  structures  from  laminated  composites  are  supposed  to  play  an  increasing  role  in  future, 
especially  for  aerospace  applications.  The  thermal  loads  will  in  many  cases  decisively  influence  the 
design.  For  example,  this  applies  to  the  load-carrying  structure  of  spaceplanes,  large  orbital  struc¬ 
tures  or  even  primary  components  of  tomorrow's  commercial  airliners. 


2  THERMAL  ANALYSIS  WITH  MSC/NASTRAN,  VERSION  66A 

Subsequently,  the  most  important  capabilities  required  from  a  Finite-Element  thermal  analysis  code 
are  specified  and  compared  with  the  features  of  MSC/NASTRAN’s  version  66A  [1],  Some  of  the 
following  topics  are  specific  for  composites,  others  refer  to  thermal  analysis  in  geni  .1. 

•  Suitable  elements  taking  into  account  the  thermal  anisotropy 
and  inhomogeneity  of  laminated  composites 

MSC/NASTRAN  provides  2D-  and  3D-elemems  with  thermal  anisotropy.  However,  there 
is  no  element  for  laminates.  This  shortcoming  has  been  tackled  by  our  elements  QUADTL 
and  QUADQL,  which  are  described  in  chapter  3  of  this  paper. 

•  Boundary  conditions  and  loads 

They  should  be  as  flexible  as  possible,  e.g.  they  should  include  time  and  temoeraiure 
dependent  convection,  temperature  dependent  emissivity,  etc. 

MSC/NASTRAN  has  full  temperature  but  no  time  dependent  convection.  The  emissivity 
can  only  be  made  temperature  dependent  in  transient  analysis,  not  in  steady  state  calcu¬ 
lations. 

•  Nonlinearities 

They  occur  in  most  thermal  problems  duo  to  radiative  boundary  conditions  and  tempera¬ 
ture  dependent  thermophysical  propterties.  Especially,  the  latter  is  important  for  compo¬ 
sites,  since  the  conductivity  in  fibre  direction  of  an  ordinary  carbon-epoxy  increases 
by  about  60%  between  30  "C  and  120'C  .  while  alumina  (99%)  has  a  nearly  constant 
conductivity  in  the  range  of  O'C  to  500  “C  (2], 

MSC/NASTRAN  can  handle  radiative  boundary  conditions,  while  temperature  dependent 
conductivities  in  transient  analysis  are  difficult  to  incorporate  and  temperature  depen¬ 
dent  heat  capacities  are  Impossible. 


View  factors  for  radiation  exchange 


This  is  an  important  issue  for  many  aerospace  structures,  especially  satellites. 
MSC/NASTRAN  has  a  module  for  automatic  calculation  of  these  factors.  However,  it 
works  very  unsatisfaciory  at  present. 

•  User  friendliness 

Especially  modelling  of  boundary  conditions  should  be  easy. 

MSC/NASTRAN  has  special  heat  boundary  elements  for  convective  and  radiative  bound¬ 
ary  conditions.  They  are  supported  by  the  interface  program  PATNAS  (inierface  to  PA- 
TRAN)  [3] .  However,  modelling  of  temperature  boundary  conditions  in  trar  sient  analysis 
costs  a  lot  of  effort  .  Another  inconvenience  is  the  lack  of  an  automatic  time  stepping 
algorithm  in  transient  problems. 

•  Interfaces  to  pre-  and  postprocessors 

Items  like  heat  flux  loads  and  radiative  as  well  as  convective  boundary  conditions  were 
very  poorely  supported  up  to  version  2.2  of  PATNAS.  This  has  been  improved  recently 
within  version  3.0  (January  1991)  (31.14). 

According  to  MacNeal-Schwendler  the  shortcomings  mentioned  above  will  be  eliminated  in  version 
68,  which  is  supposed  to  be  released  in  1993. 


3  THERMAL  ANALYSIS  OF  LAMINATED  COMPOSITES 
VIA  THE  FINITE  ELEMENTS  QUADTL  AND  QUAOQL 

The  thermal  analysis  of  composites  becomes  an  anisotropic  and  inhomogeneous  problem,  if  fibres 
and  matrix  exhibit  different  thermal  conductivities.  Especially,  this  applies  to  high-  modulus  carbon 
fibre  reinforced  plastics,  where  the  degree  of  thermal  anisotropy  can  reach  180.  This  is  the  ratio 
of  the  principal  conductivities  of  a  lamina  in  fibre  direction  and  perpendicular  to  it  ( k|/ky ) . 

A  Thermal  Lamination  Theory  (TLT)  applying  to  plate  and  shell  structures  with  high  degree  of  thermal 
anisotropy  has  been  developed  (5], (6].  It  assumes  a  linear  temperature  distribution  over  the  thick¬ 
ness  and  can  be  regarded  as  the  thermal  analogue  to  the  Classical  Lamination  Theory  (CLT) .  The 
stacking  sequence  of  the  laminate  is  taken  into  account  by  introducing  'moments  of  heat  flux*  as 
equivalent  forces.  The  isoparametric  quadrilateral  finite  element  QUADTL  based  on  the  TLT  has  eight 
degrees  of  freedom  -  one  at  each  node  for  the  temperature  of  the  laminate's  reference  surface 
and  another  for  the  temperature  gradient  in  thickness  direction.  Steady  state  calculations  of  plates, 
subjected  to  a  very  concentrated  heat  flux,  showed  a  favourable  agreement  with  fully  three-dimen¬ 
sional  calculations,  while  the  computing  time  was  reduced  by  two  orders  of  magnitude.  However, 
the  heat  flux  calculated  from  the  derivatives  of  the  temperature  field  is  approximated  worse.  More¬ 
over  the  linear  temperature  distribution  over  the  thickness  might  be  to  rough  in  the  case  of  transient 
heat  conduction. 

That's  why  a  higher  order  theory,  called  the  Quadratic  Lamination  Theory,  has  been  established 
[7] .  The  finite  element  QUAOQL  is  related  to  the  new  theory  and  has  twelve  degrees  of  freedom. 
Numerical  examples  provided  excellent  temperature  and  heat  flux  results  for  steady  state  as  well 
as  transient  problems. 
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simulation  of  the  Thermoforming  Process  of  Continuous  Fibre- 
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R.  Scherer,  Technical  University  Hamburg-Harburg ,  Germany 

RESEARCH  PROJECT 

Thermoforming  is  one  of  the  most  promising  techniques  to 
manufacture  components  made  of  continuous  fibre  reinforced 
thermoplastic  composites.  The  main  feature  is  the  molten  state  of 
the  thermoplastic  matrix  during  forming.  The  laminate  can  be 
modelled  as  an  alternating  sequence  of  fibre  rich  layers  and  resin 
Interlayers,  where  the  height  of  the  interlayer  is  approximately 
20  nm  for  a  carbon  fibre/polypropylene  composite.  The  individual 
fibre  rich  layers  have  to  slip  onto  each  other  to  accomodate  onto 
a  rigid  tool. 

"Ply-Pull-Out"  experiments  characterized  the  "interply-slip" 
behaviour  in  dependency  on  the  interlayer  height,  forming  speed, 
processing  temperature,  pressure  applied  normal  to  slip  direction 
and  lay-up  configuration.  These  experiments  represent  a 
one-dimensional  characterization  of  the  interply-slip  effect.  The 
measured  data  deal  as  input  parameters  for  a  Finite  Element 
Analyses  of  a  therrooforming  process. 

First,  the  material  law  was  defined  by  two  different  ways  with 
MARC,  and  the  "ply-pull-out"-experiments  were  simulated,  whereby 
the  material  law  is  nonlinear  and  rate  dependent. 

Two-dimensional,  8-noded  plane  strain  elements  were  used  for  the 
modelling  studies. 

This  material  nonlinear  behaviour  is  superimposed  by  a  geometrical 
nonlinear  mode  in  two-dimensional  modelling  studies.  Initially 
flat  laminates  were  thermoformed  into  a  90 '-angle  and  a  "top-hat"- 
section  in  between  two  rigid  dies.  The  deformation  behaviour  of 
two-layered  laminates  were  investigated  so  far,  and  the  resin 
interlayer  was  reduced  to  a  contact  surface  in  between  deformable 
elements.  A  user-defined  subroutine,  written  in  "Fortran",  defined 
the  properties  of  this  interlayer.  The  main  advantage  of  a 
computer  simulation  is  the  posslblity  to  vary  the  parameters  of 
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the  interlayers  so  that  the  stress  transfer  from  the  interlayers 
into  the  fibre  rich  layers  can  be  visualized.  This  was  not 
possible  during  thermoforming  experiments  which  have  been  carried 
out  simultaneously.  The  more  the  viscosity,  namely,  the  "inner 
friction"  of  the  resin  interlayers  is  increased,  the  higher  the 
tendency  towards  buckling  of  plies  is  pronounced. 

Both,  the  incremental  deformation  behaviour  during  experimental 
and  computer  simulation  studies  lead  to  the  same  results  which 
encourages  to  say  that  an  appropriate  material  model  for  the 
interlayer  was  found. 

Three-dimensional  modelling  studies  are  currently  in  process.  A 
therrooforming  process  of  a  hemispherical  part  in  between  two  rigid 
dies  is  modelled  with  shell  elements.  The  material  properties  of 
the  fibre  rich  layers  are  strongly  anisotropic.  On  the  one  hand 
side,  the  fibres  dominate  the  elastic  response  against  deflection. 
On  the  other  side,  the  matrix  possesses  a  diminishing  material 
response  in  the  direction  perpendicular  to  the  fibres. 


MARC  Is  able  to  simulate  material  and  geometrical  nonlinear 
behaviour  as  described  above.  The  main  drawbacks  result  from  a 
limited  usage  of  different  program  facilities  simultaneously.  For 
example,  the  material  nonlinear  description  requires  a  time 
incrementation  for  the  calculation  of  rate  depenedent  properties 
and  the  geometrical  nonlinear  analysis  requires  a  time 
incrementation  for  the  movement  of  the  dies.  Unfortunately,  both 
analyses  cannot  be  carried  out  at  the  same  time.  Therefore,  it  is 
not  possible  with  MARC  to  calculate  a  forming  process  of  a  rate- 
dependent  visco-plastic  material. 

A  general  disadvantage  of  some  FE-programs  is  the  "suspicious" 
definition  of  convergence  criteria  for  a  nonlinear  iterative 
calculation.  Additionally  it  is  recommendable  that  the  size  of 
increments  should  be  automatically  adapted  to  the  current  program 
in  the  ideal  case.  If  the  increments  are  manually  defined,  both 
small  and  large  increments  can  lead  to  significant  different 
results.  Of  course,  a  rough  incrementation  should  be  used  to  save 
computing  time,  but  it  often  leads  to  a  large  number  of  recycles 
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within  one  increment.  Therefore,  user  manuals  should  outline  the 
influence  of  the  different  paramerers  in  a  nonlinear  analysis  in 
more  detail. 


Abstract  submitted  for  oral  presentation  on  the  workshop  "Finite 
Element  Methods  for  Composites",  5/6  September  1991,  Imperial 
College  of  Science,  Technology  and  Medicine,  Centre  for  Composite 
Materials,  London,  England. 
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APPENDIX  B 
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ABSTRACT 

A  procedure  for  predicting  the  properties  of  composite  laminates  in  three 
dimensions  is  presented.  Analysis  has  been  carried  out  using  two 
commercially  available  finite  element  analysis  codes.  Constraint  conditions 
have  been  applied  to  the  model  to  simulate  the  conditions  of  classical  laminate 
theory  and  the  in-plane  predictions  of  both  procedures  are  compared.  Mesh 
refinement  in  both  the  thickness  and  in-plane  dimensions  has  been  investigated 
and  optimised.  The  resulting  procedure  will  predict  out-of-plane  as  well  as  in¬ 
plane  properties  of  composite  laminates  using  economically  small  models. 


I  INTRODUCTION 

The  first  use  of  advanced  continuous  fibre  reinforced  composites,  in  the 
early  1960s,  was  in  direct  substitution  for  metallic  panels  on  aircraft.'  At 
this  initial  stage  the  laminates  were  mechanically  fastened  to  a  metallic 
framework.  The  great  reduction  in  weight  was  sufficient  to  justify  the  high 
cost  of  the  materials  and  the  manually  intensive  fabrication  procedures 
employed. 
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In  this  plate-like  form  it  was  justifiable  to  assume  the  simplifications 
associated  with  the  engineering  mechanics  approach  applied  to  thin  plate 
behaviour.  Essentially  this  approach  neglects  through  thickness  behaviour 
of  the  material  including  shear  stresses.  In  addition,  linear  elastic  behaviour 
was  assumed  since  the  materials  were  based  on  thermoset-based  matrix 
systems  with  limited  strain  capacity.  Consequently  a  design  methodology 
emerged  which  is  heavily  based  on  these  assumptions.  This  allowed  use  of 
the  existing  metal-based  analysis  procedures,  after  minor  modifications  to 
allow  for  orthotropy  of  material  properties,  for  the  analysis  of  composite 
structures. 

This  methodology  relies  on  measured  in-plane  mechanical  properties  of 
the  basic  composite  materials.  For  a  unidirectional  composite  these  arc: 
two  elastic  moduli,  one  shear  modulus  and  one  lateral  contraction  ratio, 
Poisson’s  ratio.  Procedures  based  on  classical  laminate  theory  (CLT)‘  are 
employed  for  prediction  of  in-planc  elastic  properties  of  composite 
laminates  containing  plies  oriented  at  any  angle  (multi-angle  laminates). 
Effective  elastic  properties  derived  from  CtT  are  then  used  in  engineering 
mechanics  approaches  which  enable  the  analysis  of  simple  plate-like 
structures  subjected  to  selected  idealised  loading  modes.  The  elastic 
properties  derived  from  CLT  are  also  used  in  modified  numerical  analysis 
tools  such  as  the  finite  element  analysis  (FEA)  method;  this  enables  the 
analysis  of  structures  of  complex  geometry. 

In  recent  years  many  developments  and  changes  have  taken  place. 
Composite  materials  with  significant  toughness  have  been  introduced,  first 
in  thermoplastic-based  form*’*  and  later  thermoset-based.'  *  These 
materials  have  a  significant  resistance  to  delamination,  the  traditional 
weakness  of  composites.  The  ability  to  join  composites  to  other  materials 
and  to  one  another  has  improved  significantly  through  the  development  of 
better  adhesives  for  all  composites,  and  bonding  and  welding  techniques 
suitable  for  thermoplastic-based  composites.*  Fabrication  techniques 
capable  of  achieving  rapid  forming  of  complex  shapes,  such  as  rubber  block 
stamping,  have  been  developed  for  thermoplastic  composites.’  Methods 
applicable  to  thermoset-based  composites  are  also  being  developed.  Finally 
and  significantly,  experience  with  composite  materials  indicates  that  the 
majority  of  failures  which  occur  below  the  design  load  level  are  due  to 
through  thickness  stresses.' 

In  view  of  these  developments,  composite  structures  arc  no  longer 
restricted  to  simple  plate  form  attached  to  metallic  frames.  Plate-like 
laminates  can  be  formed  into  complex  shapes  with  intricate  integral 
supporting  frameworks.  Composites  are  joined  together  to  form  all- 
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composite  structures.  In  addition,  modern  composite  materials  are  used  to 
build  three-dimensional  solid  components.’  *® 

The  use  of  existing  analysis  methods  and  their  assumptions  for  present- 
day  materials  and  applications  must  be  questioned.  Analysis  procedures 
and  tools  which  account  for  through  thickness  behaviour  of  the  materials 
are  needed.  Investigators  have  proposed  procedures  based  on  three- 
dimensional  constitutive  equations"  and  the  longwave  approach.'*  We 
have  proposed  a  numerical  approach  based  on  the  finite  element  analysis 
(FEA)  method  which  allows  prediction  of  both  in-plane  and  out-of-plane 
properties  of  composite  laminates.'^  In  this  paper  we  extend  the  approach 
and  consider  the  effect  of  finite  element  mesh  density  on  the  predicted 
value'  We  compare  thercsults  to  those  produced  using  CLT,  and  propose  a 
cost-efficient  procedure  for  predicting  through  thickness  properties  as  well 
as  in-plane  properties  without  the  restrictions  of  plate  behaviour.  This 
procedure  can  be  used  for  the  investigation  of  the  effects  of  inclusion  of 
interlayers  and  laminae  of  different  thicknesses  in  a  composite,  as  well  as 
examination  of  detailed  stress  distribution  within  the  laminate. 


2  MATERIAL 

The  material  studied  is  a  multi-angle  laminate  of  continuous  carbon  fibre 
reinforced  thermoplastic  composite  composed  of  semi-crystalline  poly- 
etheretherketone  (PEEK)  thermoplastic  polymer  matrix  reinforced  with 
61%  by  volume  of  continuous  Hercules  AS4  carbon  fibres.  A  procedure  for 
determining  the  properties  of  the  basic  unidirectional  material  has  been 
presented  in  Ref  14.  Measured  properties  inherently  include  an  amount  of 
uncertainty.  An  analysis  of  the  sensitivity  of  the  present  predictive 
technique  to  the  accuracy  of  the  input  data  has  been  carried  out  and  will  be 
published  elsewhere.  For  the  present  modelling  a  consistent  set  of  input 
data,  shown  in  Table  I,  has  been  used.  The  conventional  axes  nomenclature 
is  used,  namely  direction  I  being  the  fibre  direction  and  directions  2  and  3 

Table  I 

Material  properties 

£,  fj  C,j  Cjj 

(GPa)  (GPa)  (GPa)  (GPa) 

1300  9  7  5  51  3  59 
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being  the  two  orthogonal  directions,  in-plane  and  oui-of-plane  respec¬ 
tively.  The  lamina  is  assumed  to  exhibit  transverse  isotropy,  that  is  the 
properties  in  directions  2  and  3  are  assumed  to  be  identical.  The  results 
presented  here  are  for  the  analysis  of  a  quasi-isotropic  laminate  [-h45/0/ 
-45/90],.  As  described  below,  the  method  of  analysis  used  is  suitable  for 
any  symmetric  balanced  laminate. 


3  FINITE  ELEMENT  ANALYSIS 


Finite  element  modelling  has  been  carried  out  using  two  commercially 
available  general-purpose  codes,  LUSAS  and  ABAQUS.  Most  analyses 
use  20-noded  isoparametric,  three-dimensional  solid  elements.”  Two 
analyses,  for  comparison  with  the  results  from  CLT  were  carried  out  using 
cight-noded  elements.  Orthotropic  material  properties  were  used;  all 
material  behaviour  was  assumed  to  be  linear  elastic. 

The  development  of  the  finite  element  model  is  shown  in  Fig.  I.  The 
analysis  is  for  an  infinite  laminate,  which  is  represented  by  a  block  as 
shown.  Since  the  laminate  is  balanced  and  symmetric,  it  exhibits  no 
coupling  of  overall  deformation  either  in-plane  to  out-of-plane  or  axial  to 
shear.  Since  the  laminate  is  represented  by  the  block,  the  deformation  of  the 
representative  block  is  assumed  to  be  of  the  same  form. 

The  representative  block  has  three  planes  of  symmetry,  shown  by  the 
dashed  lines  in  Fig.  I.  Thus  the  deformation  of  the  overall  block  can  be 
analysed  via  the  deformation  of  an  J  of  it,  the  cuboid  ABCDx/lytS,  The 
deformation  of  ABCDct^y^  for  the  imposition  of  in-plane  stress,  in  the  >•- 
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direction,  is  shown  in  Fig  2.  The  stress  is  applied  via  a  prescribed 
displacement  to  the  nodes  of  the  face  CDiSy.  The  planes  AB/li  and  AD^i 
are  restrained  to  remain  in  their  original  planes  from  the  symmetry  shown 
in  Fig.  I.  Constraint  equations  are  applied  to  the  nodes  on  planes  BCy/f  and 
forcing  them  to  remain  parallel  to  their  original  directions.  Thus  the 
deformation  of  a  balanced  symmetric  laminate  is  imposed  on  the 
representative  block. 

Results  from  finite  element  analysis  include  reactions  to  earth.  R,  at 
surface  nodes,  and  the  displacement  of  all  nodes.  A.  The  elastic  constants 
for  the  in-plane  loading  in  the  y-direction.  shown  in  Fig.  2,  are  calculated 
from  the  finite  element  results  via 

'  '&,/!,  A,//,  A,//, 

The  elastic  constants  from  the  other  loading  cases  are  calculated  from  the 
reactions  and  displacements  from  the  respective  finite  element  analyses  by 
analogous  calculations. 

The  finite  element  grids  varied  in  the  number  of  elements  both  for  the 
O'lt-of-plane,  ?-direction,  and  the  ir.  plane  directions,  x  and  y  (Fig.  2).  The 
shape  of  the  grid  in-plane  was  always  square,  that  is/,  =  /,(Fig.  2).  Analyses 
were  carried  out  using  cubic  elements,  and  with  longer  aspect  ratio  elements 
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Table  2 

Finite  element  griilt 


Number  of  elements 
(in*plane) 

Shape 

Lengih/lhickness 

Package  used 

1  X  1 

Cubic 

0083  33 

LUSAS  and  ABAQUS 

2x2 

Cubic 

01567 

LUSAS 

5x5 

Cubic 

04157 

LUSAS  and  ABAQUS 

5x5 

Long 

2083 

ABAQUS 

10  X  10 

Cubic 

0833  3 

ABAQUS 

lOx  10 

Long 

4- 16? 

ABAQUS 

I5x  15 

Cubic 

1  25 

ABAQUS 

15x  15 

Long 

625 

ABAQUS 

with  the  two  in-plane  lengths  being  five  times  the  out-of-plane  length,  so  the  I 
ratio  IJI,  is  increased.  This  ratio  is  designated  the  length/thickness  ratio,  l/i.  ' 

Finite  element  analyses  were  carried  out  for  various  column  heights  with  a 
single  element  in-plane;  asdesenbed  below,  these  results  allowed  deduction  . 
of  the  number  of  elements  required  per  lamina  in  the  thickness  direction. 
Further  analyses  were  carried  out,  for  the  given  column  height,  for  various  ‘ 
numbers  of  elements  in-plane  up  to  IS  x  15  elements.  The  grids  analysed 
using  20-noded  elements  and  12  elements  out-of-plane,  that  is  thickness, 
are  shown  in  Table  2. 

4  RESULTS 

4. 1  Effect  of  Through  Thickness  Density 

The  effect  of  through  thickness  density  was  investigated  by  the  analysis  ! 

of  a  column  containing  a  single  cubic  element  in-plane.  The  finite  element  ! 

package  used  was  LUSAS,  using  20-noded  elements.  The  laminate  was  I 

analysed  using  grids  of  8, 1 2  and  24  elements,  that  is  with  2, 3  and  6  elements  i 

per  lamina  lespcctivcly  As  expected,  the  results  for  the  two  in-plane 
directions,  for  values  of  both  Young’s  modulus,  E,  and  Poisson’s  ratio,  v,  ' 

were  identical.  The  results  are  summarised  in  Table  3.  These  results  show  | 

that  three  elements  per  lamina  allow  accurate  prediction  of  the  through  ' 

thickness  modulus,  f.;  the  value  is  identical  to  that  obtained  using  six  i 

elements  per  lamina  There  is  a  small  but  measurable  discrepancy  in  the 
values  of  Poisson’s  ratio,  v.^;  this  discrepancy  is  less  than  03%  and  is 
therefore  considered  insignificant.  For  the  prediction  of  in-plane  stiffness 
the  results  in  Table  3  show  that  two  elements  per  lamina  may  be  sufficient, 
although  this  could  lead  to  some  discrepancy,  around  08%,  in  the 
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TABl.t  3 

Effect  of  I  otunm  hci^hi 


In-piunc.  VA  stress 

*’a.  *'«• 

elcmcnls  per  lamina 

(GPa) 

2 

4777 

0  218  5  0  29SV 

3 

4770 

0^2167  02965 

Oul-of*plunc.  -T  stress 

E. 

elcmcnls  per  l.imina 

(GPa) 

") 

10-12 

0  062  5H 

3 

1065 

0066  IK 

6 

10  65 

0  066  37 

prcdicled  value  of  the  in-plane  Poisson's  ratio,  It  was  concluded  that  the 
present  analysis  should  proceed  using  three  elements  per  lamina.  All  results 
presented  below  used  grids  with  12  elements  in  the  thickness  direction,  the 
c-dircetion,  representing  four  composite  plies. 

4.2  Comparison  of  Results  from  Finite  FIcmcnl  Analysis  and  Classical 

Laminate  Theory 

T.'.e  ''■uislrainls  of  classical  laminate  theory  (CLT)  can  be  imposed  m 
finite  element  analysis  using  a  single  column,  that  is  one  cubic  element  in- 
plane.  of  eighi-noded  elements.  The  allowable  deformation  of  such  a 
column  of  these  elements  imposes  the  condition  of  zero  through  thick  ness 
stress,  that  is  the  condition  of  plane  stress  assumed  in  CLT 

T*hi.i  4 

Cfunpun^oii 

11/  /rmre  i 

[•U’tuvnt  onJ  Cl.T  r<'w«/r\ 

In-planc.  .vv  stress 

E. 

'  «»  * « • 

(GPa) 

CLT 

51  08 

030.1 

ABAQUS 

5108 

0302  7  0  267  1 

LUSAS 

5108 

0302  7  0  267  4 

Oul-of-planc.  ::  stress 

(GPa) 

V., 

ABAQUS 

1073 

0056  I  5 

LUSAS 

1073 

0056 16 
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The  results  for  classical  laminate  theory  were  calculated  using 
COMLAN,  an  implementation  of  classieal  laminate  theory  developed  in- 
house  (ICI)  using  the  LOTUS  123  spreadsheet  The  finite  element  results 
were  predicted  using  both  packages,  LUSAS  and  A  BAQUS.  The  results  are 
shown  in  T'bl?  4.  These  results  show  complete  agreement  between  all 
methods  of  prediction.  Classical  laminate  theory  cannot,  by  definition, 
predict  any  results  out-of-plane,  the  r-direction.  Finite  element  analysis 
using  the  same  condition,  plane  stress,  as  CLT,  predicts  identical  in-planc 
results  to  CLT.  and  also  predicts  out-of-planc  results  for  this  plane  stress 
condition. 

4.3  Results  from  Larger  Finite  Element  Grids 

Finite  element  analyses  were  earned  out  for  larger  size  grids,  using  both 
cubic  elements  and  elements  with  longer  aspect  ratios,  as  shown  in  Table  2 
The  predictions  for  the  application  of  in-plane  stress  arc  shown  in  Figs  3  5 
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All  results  are  plotted  against  the  ratio  of  the  in-plane  length  to  thickness, 
//r.  The  value  of  this  ratio  indicates  the  relative  influences  of  the  edge  effect 
on  the  result.  For  comparison  the  graphs  also  include  the  predictions 
assuming  plane  stress,  described  in  Section  4.2. 

The  results  show  excellent  agreement  for  the  two  packages  used;  the 
maximum  discrepancy  occurs  for  the  value  of  the  in-plane  Poisson's  ratio. 

predicted  using  a  single  cubic  column;  the  discrepancy  is  2%.  The 
results  from  the  different  grids  all  lie  on  smooth  curves.  These  curves  are 
discussed  in  detail  in  Section  5. 

The  application  of  out-of-plane  stress  to  the  difTercni  grids  produced 
insignificant  variation  in  the  predicted  value  of  the  predicted  value 
varied  between  10  60  and  10-67  CPa  within  the  range  of  the  length 
thickness  ratio  shown  in  Figs  3-5.  Thus  prediction  of  the  oui-of-plane 
modulus  is  independent  of  the  value  of  the  lengih/thickness  ratio  and  is  not 
significantly  difTereni  from  the  value  predicted  for  plane  stress  conditions. 
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10  73  GPa.  shosvn  in  Table  4.  The  prediction  of » is.  however,  dependent 
on  the  ratio:  the  variation  is  shown  in  Pig  6.  The  value  obtained  assuming 
plane  stress  is  also  included  for  comparison  These  results  are  similar  to  the 
results  for  »•,.  (Pig.  5).  with  all  results  on  the  same  curve.  This  trend  is 
discussed  in  detail  in  Section  5. 


5  DISCUSSION 

The  results  shown  in  Figs  3  6  may  all  be  described  in  a  single  fashion.  The 
results  from  the  different  grids  may  be  plotted  on  single  curves.  In  other 
word>.  the  value  predicted  is  related  to  the  length  thickness  ratio  of  the  grid. 
The  value  of  this  ratio  represents  the  relative  inlUience  of  the  edge  effects  on 
the  predicted  result.  Assuming  that  the  si/c  of  the  edge  effc'ci  region  is 
proportional  to  the  laminate  thickness,  then  as  the  length,  thickness  ratio 
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tncreasos  (he  influence  of  ihc  edge  region  diminishes.  Furihcrmore.  the 
material  within  the  edge  region  is  less  constrained  than  that  in  the  interior 
region  of  the  laminate,  and  hence  apparently  less  siilT.  As  described 
previously,  the  overall  distortion  of  the  grid  has  been  set  for  a  balanced 
symmetric  laminate,  allowing  no  overall  non-uniform  distortion,  as  shown 
in  Fig  2.  However,  it  is  clear  that  distortions  must  occur  within  the  laminate 
at  the  interfaces  between  laminae  oriented  in  different  directions  The 
mathematical  formulation  of  a  single  column  of  eighi-ntxied  elements 
allows  no  such  distortions  to  occur  The  amount  of  permitted  distortion 
increases  as  the  model  compliance  is  increased  until  a  value  of  /  r  is  reached 
where  the  entire  model  represents  the  edge  effect  icg.on  Beyond  this  vdluc 
of  /  f  further  increases  in  the  model  si/e  represent  increased  interior  region 
and  hence  increased  stiffness. 

The  curves  for  the  in -plane  modulus,  F.\.  and  the  in-planc  PoiNs»>n's  ralnx 
*'tr  (Figs  3  and  4),  follow  similar  irenus  The  maximum  value  is  predicted 
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uMng  the  assumplions  of  plane  stress;  the  prediction  is  a  minimum  for  a 
length, thickness  ratio  of  about  I.  and  then  increases  Both  curves  tend 
towards  a  constant  value  which  lies  between  the  value  predicted  for  plane 
Stress  and  the  value  predicted  for  the  single  column  of  cubic  20-nodcd 
elements  Thus  the  constraints  imposed  by  plane  stress,  allowing  no 
distortion  between  the  laminae,  increases  the  apparent  stiffness  and  the 
overall  lateral  contraction;  this  is  the  result  expected  intuitively.  Allowing 
iiitcr-laminae  distortions  to  occur  produces  a  rapid  reduction  in  the 
apparent  stiffness  and  lateral  contractions,  reaching  a  minimum  value  as 
the  in-planc  thickness  equals  the  thicknes ..  At  this  stage  the  distortion  is 
that  of  an  edge  region.  As  the  //r  ratio  is  increased,  the  additional 
deformation  is  that  of  a  more  constrained  interior  region,  thereby 
increasing  the  total  effeetivc  stiffness  to  a  value  between  that  predicted  for 
plane  stress  and  the  smallest  length, thick  ness  ratio,  the  single  cubic  column 
of  20-noded  elements. 

The  results  for  the  two  values  of  foisson  s  ratio  relating  displacements 
Ill-plane  and  out-of-plane,  v,.  and  v.^  (Figs  5  and  6),  show  similar  trends. 
The  predicted  values  are  minima  for  the  condition  of  plane  stress,  when  no 
internal  distortion  is  allowed.  The  lack  of  distortion  would  decrease  the 
lateral  contractions  in  one  plane  arising  from  the  application  of  stress  in  the 
other  plane.  The  predicted  values  reach  maximum  values  at  length.' 
thickness  ratios  of  about  l.and  then  decrease  with  the  predicted  value  at  an 
infinite  length 'thickness  ratio  apparently  lying  between  the  value  for  plane 
stress  and  the  value  predicted  for  the  single  column  of  20-noded  cubic 
elements  An  analogous  argument  as  used  above  explains  this  behaviour 

Conducting  the  complete  series  of  models  permits  accurate  evaluation  of 
the  mechanical  properties.  However,  the  models  of  the  high  /,  /  ratio  are 
costly  Therefore,  in  the  interest  of  economy,  it  is  useful  to  establish  two 
hounds  within  which  the  predicted  properties  will  lie  One  bound  is 
obtained  using  the  assumptions  of  classical  laminate  theory,  and  the  other 
hound  is  obtained  from  finite  clement  analysis  of  a  single  column  of  20- 
noded  cubic  elements.  However,  the  through  thickness  Young's  modulus. 
F...  may  be  obtained  accurately  from  the  result  using  the  assumptions  of 
classical  laminate  theory,  that  is  the  finite  element  analysis  of  a  single  cubic 
column  of  eight-noded  elements.  The  results  for  the  laminate  considered 
here  are  summarised  in  Table  5.  The  predicted  values  have  been  derived 
from  an  examination  of  the  curves  shown  in  Figs  3  6. 

The  bounds  shown  in  Table  5  are  generally  reasonably  narrow.  The 
exception  is  the  in-plane  Poisson's  ratio,  v,,;  the  variation  between  the 
bounds  IS  more  than  30“/.,  This  behaviour  is  not  surprising  since  this 
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mochanjcal  properly  is  expected  to  be  highly  sensMive  to  mter-lanun.ie 
di-^lornons.  A  more  accurate  value  for  this  parameter  could  onlv  be 
obtained  via  finite  element  analysis  of  a  large  grid  with  long  aspect  ratio 
elements.  Reasonable  bounds  for  all  other  parameters  may  he  predicted 
using  the  single  columns  as  shown  in  Table  5. 

The  predicted  value  of  the  in-plane  modulus.  49  GPa.  compares  well  wjih 
(he  measured  value  of  48  6GPa.‘‘‘  The  measured  value  of  the  m-planc 
Poisson  s  r.iiio.  i is  O  -T*"*  this  experimental  value  is  within  the  btuinds 
shown  in  Table  5  but  not  in  ativ^mcnt  the  predicted  value  I  his 
anomaly  is  the  subject  of  future  wairk. 


6  CONCLUDING  REMARKS 

The  three-dimensional  mechanical  properties  of  a  multi-angle  composite 
lanunate  have  been  predicted  using  a  I T.A-bascd  prvKcdure,  Bounds  lv»r 
elastic  constants  of  a  guasi-isoiropic  laminate  have  been  derive<l  using 
economically  small  finite  clement  grids.  Some  results  from  clasMc.tl 
laminate  theory  have  been  shown  to  be  usable.  This  procedure  should  be 
applicable  to  balanced  symmetric  laminates-  The  method  can  muv  be 
extended  to  the  prediction  of  the  clastic  constants  in  shear,  and  the  clastic 
constants  of  other  laminates. 
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